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SECTION 1 
INTRODUCTION 

1.1 BACKGROUND 

Flight Test Trajectory Control (FTTC) is an emerging, pilot-aiding 
technology valuable for improving flight test results and applicable to 
maneuver control and guidance in a number of contexts. This technique has 
provided the means for flying maneuvers consistently, precisely, and 
repeatably from flight to flight. Two versions of these controllers have 
been used: a closed-loop automatic system and an open-loop system providing 

manual piloting information. A closed-loop system used to collect perform- 
ance, pressures, and loads data from the highly maneuverable aircraft tech- 
nology (HiMAT) vehicle is described in [1], The application of the open- 
loop system on the NASA F- 111 transonic aircraft technology (TACT), F-15 
airframe/propulsion system interaction studies, and F-15 shuttle tiles test 
programs are given in [2]. 

Originally, the open-loop flight-test-trajectory guidance algorithms 
were developed on-line, in a piloted simulation using cut-and-try techniques 
that were not only man power intensive, but often produced less than optimal 
controllers. A closed-loop system designed using one-loop-at-a-time 
classical design approach is documented in [ 3 ]. Full-state feedback 
approach for closed-loop system design using linear quadratic synthesis is 
described in [4], Both these approaches have limitations in terms of design 
methodology and controller complexity. 

The work reported in reference [2] and more recent ongoing work have 
included detailed modeling of all the piloted aircraft subsystems for suit- 
able control law development. This research has identified the maneuver 
modeling as a significant technology needing further clarification and deve- 
lopment. The maneuver modeling aspect of trajectory guidance and control 
law design will be emphasized here in describing this developing technology. 
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In the present work, closed loop mechanization is carried out with the 
pilot in a supervisory role. The thrust here is in four areas 


1. Maneuver modeling, 

2. Application of modern linear multivariable synthesis, 

3. Techniques for the development of reference command and gain- 
scheduled perturbation controllers, and 

4. Exploratory investigation of the emerging nonlinear system design 
techniques via prelinearizing transforms. 


1.2 SUMMARY OF RESEARCH ACCOMPLISHED 


Maneuver modeling for all the given flight test trajectories were 
successfully carried out. In the present work 64 straight and level trims, 
51 level turn trims at a load factor of 2, and 30 level turn trims at a load 
factor of 4 were employed. A 3-D linear interpolation was employed. Thus a 
flight condition is characterized by altitude, Mach number and load factor. 
Linearized aerodynamic coefficients were also stored. This data in conjunc- 
tion with kinematic and some dynamic equations are then used to generate 
state and control histories for all the flight test maneuvers. 

Two multivariable synthesis techniques were used to obtain output 
feedback linear perturbation controllers viz. 

1) Constrained eigenstructure assignment (Shapiro & Andry, 1982) 

2) Minimum error excitation technique (Kosut, 1970) 

Out of these two, the constrained eigenstructure assignment technique was 
found to be the hardest one to iterate on primarily due to the current lack 
of understanding of the explicit relationship between eigenvectors and the 
desired time response. This difficulty, while not crucial in full state 
feedback design can become an extremely important element in output feedback 
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design. Due to its high sensitivity to the input eigenvectors, it was dis- 
carded after the initial research phase. A paper outlining these issues was 
presented at the 1985 American Control Conference (see Appendix A). The 
minimum error excitation suboptimal design was successful and is advocated 
as the main design approach. 

Exploratory research on nonlinear maneuver autopilot synthesis brought 
out the feasibility of generating the maneuver autopilots by employing sin- 
gular perturbation theory in conjunction with prelinearizing transforms. 

The methodology is outlined, though the controllers have not been worked 
out. A simple example problem to serve as an illustrative example is given 
in Appendix B. This technique appears to hold considerable promise and will 
be further pursued in future research. 

The designs obtained were scheduled as functions of Mach number and 
load factor and were tested on a linear simulation. The performance has 
been found satisfactory within the validity of the linear model assumptions. 
In the next phase these designs will be tested on a nonlinear simulation of 
the F— 1 5 aircraft. 


1.3 STUDY RESULTS 


Results from this work fall into four categories, (1) control design 
technique evaluations, (2) specific control analysis and design, 

(3) software developments for control law mechanization, and (4) specific 
control law validations. The various deliverables in these areas are: 

1) For control design technique evaluations 

a) A paper evaluating the eigenstructure assignment tech- 
nique [6] 

b) An assessment and extension of the nonlinear prelinearizing 
control technique [7], and demonstration on a simple 
example. (See Appendix B) 
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2) For control analysis and design 

a) Development and clarification of maneuver modeling equations 
beyond the analysis in references [it, 10]. 

b) Condensation of the F — 1 5 aircraft nonlinear characteri sti cs 
to a table of reference states and controls via nonlinear 
simulation trim values at approximately 145 conditions on 
the flight envelope. 

c) Decomposition of the 8 maneuvers over the flight envelope 
into 30 linear perturbation models (15 with fixed throttle 
and 15 with variable thrust). 

d) Solution of 30 output feedback gains which can be used with 
the 30 linear perturbation models to simulate maneuvers 
throughout the envelope. 

3) For software developments 

TM 

a) Extensions t^Integrated Systems, Inc. (ISI's) MATRIX^ 
SYSTEM_BUILD 1 package including a linear time varying 
FORTRAN block to give a generic linear time varying simula- 
tion capability (such a capability can be easily mechanized 
to model linear time-varying simulations of engine and air- 
craft dynamics, for example). 

b) Development of a three-dimensional interpolation program 
which converts table look-ups in altitude, Mach, and load 
factor to a one-dimensional table look-up with respect to 
time for a specific maneuver, 

c) Documented command files Aircraft-CAS model building and 
Maneuver Auto Pilot (MAP) design in the MATRIX language for 
the model generation, control law design, and simulation 
validation process, and 

4) Demonstration of the maneuver autopilot validations in a linear 

simulation. 


The linear control laws developed in this work are now ready for vali- 
dation on a nonlinear batch simulation. With suitable algebraic equation 
manipulation, the nonlinear control laws can also be mechanized on a 
nonlinear simulation. 


1.H REPORT ORGANIZATION 


Section 2 describes the aircraft and command augmentation system (CAS) 
models, the choice of outputs useful for feedback and the procedure for 
obtaining linear models. Section 3 specifies the flight test maneuvers 
analyzed in terms of the maneuver objective and how that objective can be 
reduced to a set of equations which constrain a required set of outputs. 
Section 4 outlines the linear techniques evaluated and used for the design 
of perturbation feedback controllers as well as the nonlinear tracking feed- 
back controller. Section 5 reviews the linear simulation mechanization and 
discusses the demonstration results for the eight required maneuvers. 
Conclusions are given in Section 6. The appendices describe in detail the 
issues in linear time-varying simulation, the evaluation of constrained 
eigenstructure assignment and the output error feedback controller designs, 
and nonlinear tracking control via prelinearizing transformations. 
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SECTION 2 

AIRCRAFT AND COMMAND AUGMENTATION SYSTEM (CAS) MODELS 


This section reviews the overall linear system model used to both 
develop and validate the linear perturbation control laws. Figure 2-1 below 
shows the design process used. 


rUGHT CONDITION 



i 

Figure 2-1. Use of Linear System Flight Test Trajectory Model 


The rest of the section describes the airframe, engine, and command 
augmentation (CAS) models and how they are integrated into the overall 
linear design and evaluation model. 
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2.1 AIRCRAFT MODEL 


As shown in Figure 2-1, the airframe model was obtained from the NASA 
Ames Dryden Flight Research Facility (ADFRF) LINEAR Program [11], yielding 
the aircraft model in the standard form 


x * Ax + Bu, 
y = Hx + Fu 


where 


T 

x : [6V, 6a, 6q, 60, 6B, 6p, 6r, 6<j>, 6h], 

y T : [ 6p , 6A n , 6q, 5q, 6p, ., 6r, 5r, 6M, 6a, 67, 6<)>, 66]. 


The first nine outputs are required in the linearized CAS model. 


2.2 ENGINE MODEL 


An engine model of the form 


0.2 

s + 0.2 


was assumed. However, since this lag can be compensated for with a simple 
lead-lag compensator, the problem of a slow engine time constant was not 
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included in the overall design, assuming exact cancellation. The actual 
residual due to mismodeling can best be addressed with the nonlinear simula- 
tion. Figure 2-2 below shows the overall mechanization including the CAS 
states described next. 

LEAP-LAC core EN; AT I ON 
FOE ENGINE DYNAMICS 


COMMANDED 

THROTTLE 

SETTING 



Figure 2-2. Linear Perturbation Flight Test Trajectory Control System 


2.3 COMMAND AUGMENTATION SYSTEM (CAS) MODEL 


The CAS model is highly nonlinear with gain schedules and multiplica- 
tive and saturation nonlinearities. This was linearized with a gross 
linearization and an "equivalent gain" for the multiplicative nonlinearities 
determined by the NASA ADFRF CONSEN program. The CONSEN program simulates 
the actual CAS for several time steps at a given flight condition. When all 
the transients have decayed, the ratio of inputs and outputs are then used 
to compute the equivalent gains. 
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The CAS system has no access to ailerons in the aircraft under consi- 
deration and hence, all the flight test maneuvers will be accomplished using 
throttle, elevator, rudder and differential tail. The maneuver autopilot 
outputs will be connected to the existing autopilot interface in the CAS, as 
shown in Figures 2-3 through 2-5. 
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2.4 LINEAR MODEL SUMMARY 


The overall linear model has nine aircraft states, one engine state, 
and twenty-one CAS states. The outputs to be used in the maneuver auto- 
pilot, for all maneuvers, are 


T 


y 


[6h, 6M, 6a, 6Y, 64 ., 66, 6p, 6q, 6r], 


Selection of these quantities for feedback are primarily dictated by the 
various flight test maneuvers to be executed. Thus, the altitude, Mach 
number and the angle of attack are to be either regulated about reference 
values or should be made to follow desired command histories. The choice of 
roll attitude feedback arises from a requirement to maintain wings level 
along symmetric flight test trajectories. In nonsymmetric flight test 
maneuvers, desired load factors can be sustained by an appropriate roll 
attitude command. For all flight test maneuvers, the angle of sideslip 
needs to be regulated about zero; this prompts the use of 6 feedback. An 
output feedback controller is envisaged in the present development. Since 
most output feedback synthesis approaches do not guarantee stability, it is 
decided to incorporate derivative feedbacks for as many output variables as 
feasible. This led to the selection of the body rates p, q, r to serve as 
derivative feedbacks for angle of attack, angle of sideslip and the roll 
attitude. The Mach number and the flight path angle serve as derivative 
feedbacks for altitude. Since most flight test maneuvers require angle of 
attack tracking, an integral feedback is incorporated in the angle of attack 
channel. While executing nonsymmetric flight test maneuvers, the roll 
attitude can sometimes be very close to 90°. At these conditions, small 
perturbations in the roll attitude can induce major changes in altitude and 
other states. Hence in order to provide a tight roll attitude control, two 
integral feedbacks are introduced for the roll attitude <f>. 
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SECTION 3 
MANEUVER MODELING 


The desired trajectory or maneuver to be modeled governs the subse- 
quent control law development in two ways. First, the maneuver prescribes 
the equations of motion of the vehicle on the reference trajectory and 
secondly the maneuver regimes determine the linear perturbation equations 
about the commanded trajectory from which the linear control law is 
developed. 

Aircraft flight test trajectories could be based on inertial reference 
(e.g., level-turn or 3~D guidance) or reference with respect to another 
vehicle or vehicles (e.g., in air-to-air combat). To place the flight test 
trajectory control design problem in a proper framework, the constraints 
which determine the equations of motion for various large classes of 
trajectories are described below. 

The subsections which follow give a general overview of maneuver 
modeling and the descriptions of the specified maneuvers analyzed in this 
work. 


3.1 MANEUVER MODELING OVERVIEW 


We divide single vehicle flight paths into those which require con- 
tinuous control along the trajectory and those that specify a final flight 
condition. In either case the flight test trajectory could be specified in 
terms of one of the following: 

1. Constraints on position components, 

2. Constraints on velocity components and altitude, 

3. Constraints on combinations of load, speed and altitude. 
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Combinations of these constraints could also be considered. The flight test 
maneuvers discussed in this report belong to the second and third 
categories. Reference [ H ] gives some early results in flight test 
trajectory modeling. 


3.1.1 Constraints on Position Components 


Examples of trajectories which involve position constraints along the 
flight path are 


1. i|-D guidance (x(t), y(t), h(t) are given functions of time). 

2. 3"D guidance (x, y, h are related to each other, e.g., fly along 
a hypothetical wire in space). Examples of 3"D guidance are 
approach to landing, terminal area flight paths and threat eva- 
sion for reconnaissance aircraft and bombers. This also includes 
straight and level flight and flights along predetermined paths. 

Examples of trajectories which specify position constraints at the final 
trajectory point are: 


1. il-D specification (arrive at a certain point, at a certain time, 
e.g., touchdown on the runway at a specified point). 

2. 3~D specification (fly-to-VOR, terrain following). 


Note that each of these trajectories requires position measurement. The 4-D 
guidance trajectory indirectly specifies velocity and acceleration compo- 
nents. Thus, specification of position components is the most comprehensive 
constraint on the trajectory. Such a rigid constraint may be unnecessary 
for most test maneuvers. 


3.1.2 Constraints on Velocity Components and Altitude 


While the horizontal position components do not, in general, affect 
aerodynamic variables, the altitude determines density and by itself affects 
dynamic pressure. Thus, it must always be considered as a possible variable 
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to be constrained. In fact, the altitude and dynamic pressure are so 
important that a majority of flight test trajectories will define the alti- 
tude profile (this includes maintaining constant altitude). 

Examples of this class of trajectories are: 

1. u(t), v(t), w(t) and h(t) [in other words, Mach number, dynamic 

pressure, £5(t) and a(t).] 6 ( t ) may be zero. 

2. Mach number, angle-of-attack and dynamic pressure (as in shuttle 
tile tests ) . 

Various other combinations of velocity components and altitudes could also 
be specified. 

Mach number, angle-of-attack and altitude constraints could also be 
desired at one point on the trajectory. For example, the zoom-and-pushover 
is a trajectory where angle-of-attack, Mach number and altitude are speci- 
fied at one point on the trajectory. 


3.1.3 Constraints on Combinations of Load, Speed and Altitude 

The trajectory specifications could involve components of loads along 
the three axes, velocity components and altitude. The typical load specifi- 
cation will consist of desired vertical acceleration. The desired value of 
the lateral acceleration is usually zero. The total speed is often speci- 
fied in lieu of the fore-and-aft acceleration. 

Many combinations of load, speed and altitude specifications are 
possible. Some examples are as follows: 

1. A constant load, constant Mach number level turn, 

2. A constant Mach number, constant altitude windup turn. 

Often, the desired flight trajectory for an aircraft depends upon the 
position and flight test trajectory of other vehicles. Typical examples are 
collision avoidance, air combat, or avoidance of air-to-air missiles. The 
specif ication is typically based on the position of a target aircraft with 
respect to the aircraft whose trajectory is being controlled. 
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The next section gives a specific objectives and their analytical 
development for each of the maneuvers for which autopilots were designed. 


3.2 MANEUVER MODELING 


To summarize, the objective of maneuver modeling is to generate a con- 
sistent set-of state and control histories to serve as commands and open 
loop controls for the maneuver autopilot using a data base consisting of 
trim conditions. Two sets of trims have been found adequate for all the 
maneuvers discussed here, viz, straight and level trims and level turn 
trims. To the extent feasible, kinematic relationships are employed to 
generate the state-control histories. Whenever this is not possible, 
linearized aerodynamics and engine models are used. Note that the following 
development is not restricted a particular aircraft. 

The commands and the open loop controls consist of: 

Commands: Altitude, Mach number, angle-of-attack , flight path angle, 

roll attitude, angle of side slip, roll-pitch-yaw body 
rates . 

Openloop Controls: Throttle, elevator, rudder and differential tail. 


| In the following, the maneuver modeling for individual flight test 

trajectories will be discussed in detail. It is important to note that all 

I 

| these maneuvers assume zero sideslip. 

i 
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3.2.1 Transient Trajectory 


The objective of this maneuver is to transfer an aircraft flying 
straight and level at a Mach-altitude pair to another Mach-altitude pair at 
a desired flight path angle. This maneuver is normally employed as the 
initial-terminal transient to other flight test maneuvers and hence the 
name . 

The simplest way to mechanize this maneuver to assume a cubic poly- 
nomial in time for the altitude. Thus, 


h = h 0 + a x t + a 2 t 2 + a 3 t 3 ; t 0 < t < t f 

from which 

. 2 

h = a, + 2a 2 t + 3a 3 t 


(3.D 


(3.2) 


The requirement for a cubic polynomial arises from the constraints that one 
wishes to place at the two ends, i.e., initial and final altitudes are 
specified along with altitude rates at the two boundaries. To simplify the 
development, constant acceleration along the flight path is assumed next, 
viz, 


V = V 0 + Vt, (3.3) 


where 


V - V 0 

V = , V- is the specified final speed and t f , the final 

I •*- 

time . 
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The flight path angle is readily computed from 


• - 1 A 

sin (-) 


As noted elsewhere, a data base consisting of straight and level trims 
at several Mach-altitude pairs are available. 

Assuming that the path angle is small, such that 

Lift^ weight 

along this maneuver, the angle of attack history can be generated by 
linearly interpolating between stored straight and level trim data at the 
Mach number - altitude given by the equations (3.1) and (3.3). Since the 
angle of attack will be close to the straight and level trim values, the 
aerodynamic surface setting at these trims can be used to generate approxi- 
mate open loop control settings. Assuming next that under these conditions, 
the actual drag is close to the trim values, the required thrust may be 
computed as follows. 

Assuming that the aircraft thrust is aligned with the vehicle 
longitudinal axis, the acceleration along the flight path for symmetric 
flight is given by 


V 


Tcosa 


trim 

m 


-D 


g sinY 


(3.4) 


In the expression (3.4), T is the thrust, Y flight path angle, o the 
angle of attack, D the aerodynamic drag, m the aircraft mass and V the 
velocity along the flight path. In order to compute the required thrust, 
the equation (3.4) can be manipulated to the form given below 


T 


(V + gsinY)m + D 


cosa 


trim 


(3.5) 
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To compute the throttle setting, a linear thrust-throttle characteristic 
will be assumed. Thus, 


T 

max 


T 

trim 

^trim 


(3.6) 


n. . is the straight and level throttle setting and T is the thrust 
trim max 

corresponding to maximum throttle setting. 


. . ’"'actual 


T/T 

max 


(3.7) 


During this maneuver, one expects the body rates to be small. Consequently, 
the commanded values of these quantities are zero. 

If the throttle setting emerging from this analysis is greater than 
the maximum or less than the minimum, it indicates that the assumed time of 
flight for the maneuver is unrealistic or that the model is inadequate or 
both. This quantity has then to be changed appropriately to make the 
maneuver feasible. 


3.2.2 Level Acceleration/Deceleration 

This is a wings level, constant altitude maneuver with Mach number 
constant or changing at a specified rate. The maneuver modeling for this 
trajectory is essentially a subset of the trajectory 3.2.1. Putting 
ai=a 2 =a 3 -0 in (3.1) results in 


h « h 0 , constant 

y = 0, constant (3.8) 

and V = V 0 + Vt, with V specified. 
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ALTITUDE 


The required thrust and the corresponding throttle can be computed as in 
3.2.1. The open loop control surface settings and the other commands are 
linearly interpolated from the straight and level trim data base. 

3.2.3 Pushover, Pullup 

This is a wings level, constant Mach number maneuver in which the 
angle of attack is varied a specified increment about the trim value at some 
specified rate. The maneuver is shown in Fig. 3.1. The corresponding angle 
of attack history is given in Fig. 3.2. 



Figure 3.1. Pushover - Pullup Flight Test Trajectory. 
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Figure 3.2. Angle of Attack History Along the Pushover - Pullup 
Flight Test Trajectory. 


The maneuver modeling for this flight test trajectory uses the assump- 
tion that pitch rate is close to zero. Since the aircraft is in symmetric 
flight, the flight path angle Y can be calculated as the difference between 
the pitch attitude and 0 and the angle of attack a. Thus 


Y - 0 - a, 


( 3 . 9 ) 


and a, is specified as a function of time, the flight path angle can be com- 
puted. Now, 


h - VsinY * MC SinY 


( 3 . 10 ) 
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where C is the speed of sound specified as a function of altitude h and M 

is the desired Mach number. Equation (3.10) can be analytically or 

numerically integrated to yield the altitude history. The throttle setting 

may be computed from the following equations. 

Since Mach number is to remain constant throughout the maneuver, one 

V 

can differentiate the expression for Mach number (M = — ) with respect to 
time and equate to zero to obtain 


V 

C 


V 3c; 
C 2 3h n 


0 


substituting for h from (3.10), one has 


• V 2 3c 
V " ~ 3h sinY 


(3.11) 


( 3 . 12 ) 


Equating expressions (3.12) and (3.*0, one has 


11 

C 


3c . 

•5H sln7 


TCosot ~ D 
m 


gsinY 


(3.13) 


From which 


T = [— + g] E!Sinl + 2 

L C 3h KJ Cosa Cosa 


(3.1*0 
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The drag at the commanded angle of attack can be computed from linearized 
drag coefficient specified as a function of a, the dynamic pressure and the 
reference area. As before, the linear throttle assumption is invoked to 
compute the throttle setting. 


^actual 


T . -tel"! 

T 

trim 


(3.15) 


3 . 2 . -4 Zoom and Pushover 


The zoom and pushover is a wings-level, thrust stabilized less than 1g 
maneuver. The flight trajectory is a parabolic path with the target 
Mach/altitude/angle of attack point at the apex. An illustration of the 
zoom and pushover flight test trajectory is given in Fig. 3.3. The maneuver 
begins at 0 with straight and level flight conditions. A transient-maneuver 
is performed to transfer the aircraft to the point A with all controls 
active. At the point A, the throttle is fixed at a predetermined value and 
the aircraft executes the zoom and pushover trajectory. At the point B, the 
thrust control is reinstated and a transient trajectory transfers the 
aircraft back to straight and level flight conditions at point C. 



(♦-Variable throttle Fixed throttle ^ — Variable throttle- 

Figure 3.3. Zoom and Pushover Flight Test Trajectory. 
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This is perhaps the most complicated of all the symmetric flight test 
maneuvers. This flight test trajectory consists of three phases. In the 
first phase, the aircraft is transferred from its straight and level flight 
condition to the beginning of the parabolic flight path. The second phase 
is the required zoom and pushover maneuver followed by the third phase which 
brings the aircraft back to its original straight and level flight condi- 
tions. The first and third phase maneuvers are essentially transients and 
the maneuver modeling discussed in 3.2.1 is directly applicable. The second 
phase will be analyzed in this section. 

A parabolic flight path has the following properties 

1. Horizontal acceleration is zero 

2. Vertical acceleration is constant 

3. Total energy is constant. 

Since the apex speed is specified, say V T , one has 
x « VcosY - constant = V T 

Here, x is the down range. Note that since the aircraft is in 
symmetric flight, the cross range is zero. Thus, 

-1 V T 

Y * + cos [^-] (3.16) 


a positive or negative sign has to be chosen based on whether the aircraft 
is flying towards the apex or away from it. 

From the given apex speed, angle of attack and altitude, the lift and 
drag can be computed using the straight and level trims data base using 
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linearized aerodynamic coefficients. From constant energy property, one 
has, at the apex of the parabola, 


TCosa = D 

T 


(3.17) 


From which, the throttle setting at the apex can be computed as follows 


^actual 


T . -trl 1 ? 

T trim 


(3.18) 


with lift and thrust, the vertical acceleration at the apex can be computed. 
Thus 


Tsina T + L 
m 


g = g 


a 


(3.19) 


Note that g should be a negative quantity, numerically less than the 
a 

acceleration due to gravity g. The acceleration g has to remain constant 

a 

through the parabolic path. To summarize, the aircraft trajectory approxi- 
mates the path of a projectile in a uniform conservative force field. The 
total energy of the aircraft in this field is given by 

V 2 

E = h + — - = constant (3.20) 


Expression (3.20) can be used to compute the speed along the parabolic path, 
given the altitude. 
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Next, given the altitude at which the parabola is to begin, and 
perhaps end, one can write 


h “ h 0 + h 0 t + —V 


with h 0 computed using the following relations, 


( 3 - 21 ) 


V o - /-2g a (E-h 0 ) 
-1 

Y 0 = COS 

* r\ 


V„sinY„ 


( 3 . 22 ) 


The time of flight on this parabola is easily computed with equation (3*21) 
in the general case or with the following equation if the zoom and pushover 
parabola is symmetric about its axis. 


tf - /(V|-V’) (3.23) 

Note that the throttle is to remain fixed at the value given by the equation 
(3.18). The angle of attack throughout the parabolic path is computed from 


T*a + L. + L a 


m 


- gcosY - g. 


(3.2*0 
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with T = T • n . , and small angle approximation for a has been used, 

max actual 

The open loop control surface settings are assumed to be the inter- 
polated values from the straight and level trim data base. 

3.2.5 Excess Thrust Windup Turn 


This is a maneuver with angle of attack linearly increasing from the 
wings-level trim condition to some specified final value at a specified 
rate. The maneuver is performed at constant altitude and constant Mach 
number. A schematic figure of this maneuver is shown in Fig. 3.H. 
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Figure 3.H. Excess Thrust Windup Turn Flight Test Trajectory. 


As mentioned elsewhere in this report, the trim data base consists of 
straight and level trims at several Mach-altitude pairs along with level 
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turns at several Mach-altitude-load factor conditions. Since excess thrust 
windup turn can be considered to consist of several level turn trim condi- 
tions appropriately pieced together, this maneuver model is merely a 3 - 
Dimensional interpolation using the trim data base. 


3.2.6 Constant Throttle Windup Turn 

This is a maneuver with angle of attack increasing linearly at a 
specified rate from trim to some specified final value. The maneuver is 
performed at a predetermined, constant thrust level. Mach number is main- 
tained by trading potential for kinetic energy via an appropriate altitude 
rate. 


Since Mach number is constant, one can write as in equation (3.12) 


V 


V 2 3C . v 
“ 9h SinY 


(3.25) 


Further, since the altitude, Mach number and angle of attack at the initial 
point are known, one can compute 


T - T. . /tl. . 
max trim trim 


from the trim data base. Let the throttle be fixed at a value n D . 

R 

the actual thrust 


(3.26) 


Thus , 


T « T .hr. 

R max R 


(3.27) 
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If the actual thrust T is greater than that required for level turn at the 

n 

given Mach-altitude-angle of attack condition, the Mach number can be main- 
tained constant only by a positive altitude rate. The reverse applies when- 
ever is less than T, . . Let the excess thrust over the level turn trim 

Ft trim 

be 


AT 


T - T 
R trim 


( 3 . 28 ) 


Now, one has 


Yl 

C 


3c . v 

sH 31nV 


AT » Cosa 
m 


gsinY 


(3.29) 


from which, 


y 


sin 


ATCosa 




(3.30) 


Also , 


h = VsinY 


(3.31 ) 


The expression (3.30 can be numerically integrated over one step to obtain 
the new altitude. The calculations may be repeated as many times as one 
wishes to obtain the trajectory. 
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It is important to begin this trajectory at a high-g turn, since 
fixing the throttle at a straight and level condition can result in an 
initial acceleration or a high path angle cl imb/ descent . Hence, in order to 
avoid confusion, this maneuver requires an initial and terminal maneuver. 
Thus three phases are required to execute this maneuver. 


1. A trajectory beginning at straight and level flight at the 
desired Mach-altitude pair and ending at a high-g level turn with 
all control surfaces and throttle active. 

2. Constant thrust windup trajectory. 

3. Terminal maneuver at constant altitude transferring the vehicle 
from the level turn at constant thrust conditions to straight and 
level flight. 

A typical constant thrust windup turn trajectory is given in Fig. 3.5. 
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Figure 3.5. Constant Thrust Windup Turn Trajectory Descending Flight. 


3.2.7 Constant Dynamic Pressure and Constant Load Factor Trajectory 


This maneuver is initiated at a predetermined load factor, Mach 
number, and dynamic pressure. Thus, the initiation of this maneuver is not 
necessarily the wings-level condition. This maneuver can be either an 
ascending or descending at a specified Mach number rate. Dynamic pressure 
and load factor are held constant throughout the maneuver. Altitude is 
gained or lost to maintain dynamic pressure with changing Mach number. 

As in the maneuver 3-2.6, this trajectory also requires three phases. 
The first is a level turn to achieve the required load factor from straight 
and level flight at the desired Mach-altitude pair. The third phase 
restores the aircraft to straight and level flight at the final Mach- 
altitude condition. Both these maneuvers can be constructed by 
interpolating between the stored trim data points at the current Mach- 
altitude-load factor conditions. In the following we shall discuss the 
constant dynamic pressure-constant load factor trajectory modeling. Note 
that this development does not depend explicitly on the load factor. Thus, 
this model is valid for all load factors including a constant unity load 
factor, wings level, constant dynamic pressure trajectory. 

The dynamic pressure, Q is given by 


Q = 



. v 2 


(3-32) 


Differentiating the expression (3-32) with respect to time and using the 
altitude rate equation, with Q ■= 0, one has 




p(h)v 


0 
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or 


V 


1 Ip 

2p 3h 


V 2 sinY 


(3.33) 


Since the Mach number is given by 


M = V/C(h) , 


in order to maintain a desired Mach rate, one must have 


M = 


V 

C 


V 2 

c 7 


Ih SInT 


from, which 


. v 2 9c 

v - MC + - ■& 


<3.3*0 


Equating the expressions (3*33) and (3.3^) > one has 


Y 


Sin 1 


- M C 


(1 l£ 

K 2p 3h 




(3.35) 
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There are two important conclusions that can be drawn from the expression 
(3.35) 


1. Since and are negative (density and speed of sound decrease 

with increasing altitude), a positive M will result in a climbing 
trajectory while a negative M will yield a descending path. 


2 . 


3p 3c 

This maneuver cannot be flown at altitudes where and are 
nearly zero unless the desired Mach rate is also zero. 


As before, the altitude rate equation 


h = VsinY 


(3.36) 


may be numerically integrated over a small time step to obtain the new 
altitude. The actual throttle setting is again computed from the level turn 
trim thrust and throttle setting at the current Mach-altitude-load factor as 
follows . 


T 

max 


T . /it . 
trim trim 


T = 


[g 


- — -l^v 2 ] 
2p 3h J 


msinY 

Cosa 


D_ 

Cosa 


^actual 


T/T 


max 


(3.37) 


(3.38) 


(3.39) 
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The angle of attack and drag in the expression (3*38) are the interpolated 
values from the trim data base at the current Mach-altitude-load factor. 
The Mach number as a function of time is obtained from 


M = M 0 + Mt (3. HO) 

The command body rates and open loop control surface deflections are again 
the interpolated values from the trim data base. 


3.2.8 Constant Reynolds Number and Constant Load Factor Trajectory 

This maneuver is initiated at a predetermined load factor, Mach number 
and Reynolds number. Thus, the initiation of this maneuver is not 
necessarily the wings-level condition. This maneuver can be either 
ascending or descending at a specified Mach number rate. Reynolds number 
and load factor are held constant throughout the maneuver. Altitude is 
gained or lost to maintain Reynolds number with changing Mach number. 

This maneuver model is different from 3.2.7 only in the way that one 
computes the flight path angle and the throttle setting. Hence only these 
two aspects will be discussed in the following. As in maneuver 3.2.7, the 
modeling does not depend explicitly on the load factor. Consequently, the 
following development is valid for all load factors including a wings level, 
unity load factor - Constant Reynold's number trajectory. 

The Reynold’s number, Re is given by 


VDp(h) 
y (h) 


where y is the Viscosity of atmosphere 


(3. HI) 
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Differentiating the expression (3.41) with respect to time and using the 
altitude rate equation, with Re = 0; one has 


* ■ ( i IS ■ p lh >v ‘ slny 


( 3 . 42 ) 


Equating the expression (3*42) to equation (3.34), 


y « sin 


MC 

(1 1h - 1 Ip _ 1 is-iv 2 

v y 9h p 3h C 9h ; 


( 3 . 43 ) 


As in the constant dynamic pressure, constant load factor flight test tra- 
jectory, we note that this maneuver cannot be flown at altitudes where -It, 
and are nearly zero unless the desired Mach rate is also zero. The 
actual throttle setting can be computed from the level turn trim thrust and 
throttle setting at the current Mach- altitude-load factor as follows. 


T 


r(l - 1 i£)v 2 + e i g s A ?? + —L 

y 9h p 9h' K Cosa Cosa 


(3.44) 


. T ^trim 

n actual ~ a . 

trim 


( 3 . 45 ) 


To facilitate easier computations, a FORTRAN program has been written 
to generate the required commands and open loop control settings given the 
appropriate data. A listing of this code is given in Appendix C. 
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SECTION U 

MANEUVER AUTOPILOT DESIGN 


The previous section described the required maneuver modeling, whereby 
for eight chosen maneuvers, a subset of the outputs are constrained to 
prespecified time histories. For the control analysis and design done in 
this study to have any value when applied to the F-15 nonlinear simulation 
or the actual aircraft, consistent values for all of the dynamic states and 
corresponding control values along the trajectory must be found. There are 
at least two straightforward ways to generate the required reference states 
and controls. First, one could iteratively simulate with the nonlinear 
model until an open loop law approximates the desired output time histories. 
This could be done systematically with numerical nonlinear optimization, 
using a parameterization of the control surface and thrust time histories. 

Or secondly, one can trim the nonlinear simulation at a number of conditions 
close to the desired trajectory and approximate the dynamic reference tra- 
jectory from these trim values. The latter approach was chosen both because 
of the completeness and flexibility of this tabular representation of the 
nonlinear aircraft characteristics, and because of the connection with the 
perturbation trim controllers designed to work along with the reference tra- 
jectory commands. It should also be pointed out that this "static" approach 
eliminated the need for any nonlinear simulation during the control design 
stage, apart from execution of a linearizing program [11] which contains the 
nonlinear F-15 aircraft equations. 

The next subsection shows how a table of trim values can be used to 
develop a "linear model" of the entire F-15 flight test system for design 
and evaluation of the aircraft dynamic response in specified nonlinear 
maneuvers. A second subsection outlines two different linear control design 
techniques, evaluating their strengths and weaknesses, giving perturbation 
controller designs using these two techniques. A final subsection presents 
an assessment and extension of the nonlinear prelinearizing control 
technique [12], the application of which will be demonstrated in the next 
study phase, along with linear gain scheduled controllers, on the full 
nonlinear dynamic simulation. 
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H.1 THE F-15 FLIGHT TEST SYSTEM "LINEAR" MODEL 


It must be emphasized that while the control design approaches demons- 
trated in this work are based on linear models, the maneuvers desired are 
highly nonlinear. As discussed above, the nonlinear characteristics of the 
F-15 are condensed into a table of trim values, which when properly used 
with linear perturbation controllers gives a linear time varying simulation 
which accurately represents the nonlinear aircraft response through the 
nonlinear maneuvers. Adequate control through the nonlinear maneuvers 
requires not only consistent open loop reference commands but linear 
perturbation gain-scheduled controllers as well. The gain-scheduled 
perturbation controllers are designed efficiently by decomposing the desired 
eight maneuvers into 15 linear design points for both the fixed and variable 
throttle cases. 


ll.I.I The F-15 Nonlinear Tabular Model 


The nonlinear dynamics of the F-15 can be represented over the 
envelope with a sufficient number of trim values X and U: 1 1<5 points were 
stored, distributed as shown on the altitude-Mach plane in Figure 4— 1 . 

A much coarser grid of flight conditions, than shown in this figure 
for the reference trajectory, was used for the linear perturbation models 
about the trajectory, and consequently in the linear control design stage. 
Table 4—1 shows the coarse discretization considered. Since it becomes 
increasingly difficult to trim the aircraft, particularly at high load 
factors, in the off-diagonal points in the altitude-Mach plane, only five 
conditions , the diagonal ones indicated in Table H-1 , were used with the 
three different load factors — 1g, 3g, and Hg. The eight maneuvers were 
initiated at different conditions (see Section 3) to exercise controllers at 
various points on the flight envelope. 
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Figure *4-1. F-15 Nonlinear Tabular Model Trim Point c 


TABLE *4-1. F-15 LINEAR PERTURBATION MODEL AND DESIGN CONDITIONS 


h/M 0.6 0.8 1.0 1.2 1 . *4 1.6 1 .8 2.0 


50k 
*4 Ok 
30k 
20k 
10k 


X 


X 


X 


X 


X 
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The trim conditions used in this work are those with x = 0, namely the 
net forces and moments are zero. Since altitude rate is zero, a load factor 
greater than one automatically pulls the aircraft into a level turn trim. 

The static nature of the trim points is overcome in constructing a 
dynamic reference trajectory by assuming a ramp in velocity for a specified 
time of transition between two trim points. This constant dV/dt yields the 
thrust adjustment necessary to add to the trimmed thrust for a "dynamic" 
reference thrust and V command. Conceptually, one could solve for an 
adjusted x command in all the states; however, the required computation and 
storage is large. Therefore, only the feedforward thrust (or throttle) com- 
mand is adjusted as described here, and the linear perturbation controller 
generates the extra transient control commands necessary for all the states 
to transition between the trim points. This mechanization is discussed more 
fully in the next subsection. 

In summary, the nonlinear F - 1 5 model required for eight nonlinear 
maneuvers has been represented by a tabular description of trim values of 
the nonlinear simulation at discrete points on the flight envelope. The 
number of points used to represent the reference command is much finer than 
for the linear models because the reference command effectively contains the 
nonlinear behavior in it. Load factors greater than one, level turn trims, 
have been used to generate the asymmetric models. 

By keeping the linear perturbation model grid coarse, only a few per- 
turbation controllers need be designed, reducing both linear model 
scheduling and gain scheduling requirements in the linear time-varying 
simulation validations. Nonlinear simulations in the next research phase 
will confirm whether or not the number of design points is sufficient. 
Accurately speaking, the "controller" , is not merely the linear perturbation 
gains, but the way they are mechanized along with the reference commands. 
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H.1.2 Linear Perturbation Controllers 


Currently, control systems for nonlinear plants are synthesized using 
perturbation models or the so called linearized plant models in conjunction 
with the powerful linear system design approaches. The controllers so 
obtained may be termed Linear Perturbation Controllers to denote the linear 
nature of the controllers and to indicate their function, viz, controlling 
perturbations about the reference condition. If the system is required to 
track a given command, the perturbation models need to be generated at 
several operating points along the command history and controllers designed. 
In highly nonlinear systems such as aircraft, these controllers can display 
significant variations, often requiring these to be scheduled as a function 
of the independent variable. 

The objective of the controller is to ensure that a given nonlinear 

system 


X = f(X,U) 


(4.1) 


follows a given trajectory X(t) . Here X e ! n , U e ! m . To obtain the per- 
turbation models or the linearized models, some points along the desired 
trajectory X(t ) are selected, say X^ , X^, X^ . . . . A set of controls 
corresponding to these points, , U^, • • • • are next computed such 

that 


f (X, U) = 0 


(4.2) 


Note that this is not the only possible approach. If the desired trajectory 
satisfies X * f(X,U) for nonzero X, then it can be used in the subsequent 
development. 
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To derive the perturbation model with state perturbations 6X and con- 
trol perturbations 6U, let 


X = X + 6X 
U - U + 6U 


(4.3) 


Expanding the nonlinear system (4.1) about X, U and retaining only the first 
order terms (this implies that the perturbation are small), one has 


5X = f x 5X + fySU 


(4.4) 


The subscripts denote partial derivative matrices. Note that and 
depend on X, U. The expression (4.4) describes a linear dynamic system for 
which a controller of the form 


<5U - K6X 


(4.5) 


can be designed at the operating points X^ , X 2> ... 

Next, to ensure that the system transits through these operating 
points, the following procedure is setup. 

Assume that a linear interpolation scheme between X 1 , X 2 , ' * * -n 
adequately describes the desired trajectory. Further, assume that at any 
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interpolated reference flight conditions X., in between , X ; 

f(X., u.) 2 0. In this case, the perturbation model (4.4) is given by 
J J 


6X = f 6X + f 6U - X(t) 


(4.6) 


Note that X is a piecewise constant function and appears as a disturbance in 
the perturbation model. 

As noted earlier, the perturbation controller is designed with X = 0. 
Recalling that 


X = 6X + X 


U = SU + U 


(4.7) 


X = 6X + X 


Using (4.7) in (4.6), the linearized equations can be put in the form 


X = f x (X-X) + fy(U-U) 


or 


X * f x X + f y U - Z(t) 


(4.8) 


where 
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Z(t) - f x X + fyU 



Similarly, the perturbation controller (4.5) becomes 


U - U - K(X-X) 


or 


U = K(X-X) + U 


(4.9) 


Expression (4.9) indicates the implementation of the linear perturbation 
controller, given in Figure 4-2. for clarity. 


Disturbances 



UNERR PERTURBATION CONTROLLER 


Figure 4-2. Implementation of the Linear Perturbation Controller 


Note that the gain matrix K varies as a function of some scheduled 
variable. Expression (4.8) has an interesting interpretation, viz, when the 
perturbations about the desired path are small, it describes the nonlinear 
system dynamics with a high degree of fidelity. Hence, under closed loop 
control, the dynamics given by (4.8) will be very close to that of the 
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original nonlinear system (4.1). This fact means that linear, time varying 
simulations can be mechanized in a straightforwar d way from the relation- 
ships and models discussed here. While section five reviews the concept of 
linear time-varying simulation before giving the validation results, it is 
convenient to continue our line of thought and go over the mechanization 
here in the maneuver autopilot design section. 


4,1.3 Linear Time Varying Simulation 

As noted elsewhere, for highly nonlinear systems such as aircraft, the . 
gain matrix K in (4.9) would display large variations as a function of the 
flight condition. And hence, some type of scheduling strategy will be 
essential for the satisfactory operation of the control system. To evolve 
the scheduling strategy, it is desirable to have a simulation of the system 
which has lesser complexity than the original nonlinear system. 

Examining expression (4.8) in view of the above, one finds that the 
partial derivative matrices f^ and f^ as functions of X, U have already been 
computed at the controller design stage. The reference trajectory X, U is 
also known. Since this expression describes the nonlinear system dynamics 
adequately for small perturbations , it may be used to develop a linear time 
varying simulation to evaluate the controller scheduling. Figure 4-3 gives 
the formal structure of the linear time varying simulation. This block 
diagram is structurally similar to the linear perturbation controller 
implementation given in Figure 4.2. The essential difference between them 
is that the nonlinear aircraft model has been replaced by a linear time- 
varying model with disturbance inputs. As mentioned in Section 4.1.2, X and 
U are linearly interpolated between time points. In order to simplify the 
mechanization, the partial derivative matrices f^ and f^ are stored at the 
same time points and linearly interpolated. Thus, along a contemplated 
maneuver, three or four X, U, points are chosen and the corresponding f^, 
fy, matrices are stored. Note that in the present analysis the X, U, points 
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are chosen such that at an intermediate point f(X., U.) * 0. In the course 

J J 

of simulations, if it turns out that in order to track the required X 
history, the total control required is greater than that available, it is 
indicative that either the assumed maneuver time is unrealistic or that the 
model is inadequate or both. 



Figure 4-3. Linear Time Varying Simulation 


A description of how a generic time varying simulation can be imple- 
TM TM 

mented in the MATRIX SYSTEM_BUILD model building and simulation program 
is given in Appendix D. 






H.2 LINEAR DESIGN TECHNIQUES 


Given a small set of linear trim conditions, linear perturbation con- 
trollers can be designed to track desired reference commands. The simplest 
practical controllers, once the outputs are augmented with appropriate 
integrators for desirable steady state tracking, are output feedback 
controllers. Two different output feedback approaches were evaluated, viz, 
ei genstructure assignment and the minimum error excitation technique. 


■4.2.1 Eigenstructure Assignment Design 

As described in Section 2 , the model through which the desired 
maneuver autopilots (MAP) must work includes the fully augmented aircraft: 
This actually makes the control design task more difficult than dealing only 
with the bare airframe. The combination of CAS states and integral error 
states give a considerable number of eigenvalues that are slow, and there- 
fore must be moved, and also exhibit a reasonable degree of state coupling. 
This coupling results in some sensitivity of results to the choice of 
eigenspace requested. The model used and specific results are given in 
Appendix A. Complete details of this eigenassignment approach can be found 
in [ 13 ]. It is perhaps of interest to note that in addition to the 
capability to handle controller structure constraints, this technique can be 
modified to accept partial specification of eigenvectors. This feature can 
be valuable in high order systems. In summary, to employ this synthesis 
approach, the following are required. 

(i) Based on the practical aspects of the problem, choose a minimal 
set of measurements which will permit the designer to achieve the 
desired performance. Introduce dynamic compensators such as 
integral feedbacks, lead-lag networks, etc., based on experience. 

(ii) Choose a set of desired eigenvalues and eigenvectors equal to the 
number of outputs. 
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While the selection of desired eigenvalues is often apparent in a 
given problem, the desired eigenvectors are difficult to select. Three 
approaches were developed to help guide this choice, viz, minimally restruc- 
tured eigenassignment , decoupling eigenassignment and dominant mode 
eigenassignment . 


it. 2. 1.1 Minimally Restructured Eigenassignment 

Since it is known that the closed loop eigenvectors lie in a subspace 
spanned by the columns of (AjI-A) B \ i ■ 1, ...n, linear combinations of 
these vectors were used as desired eigenvectors. The weights to be used in 
generating these linear combinations were constructed from the additional 
information that unassigned eigenvectors should be close to their open loop 
values in a least square sense. 

i».2.1.2 Decoupling Eigenassignment 

Since we are interested in having the least cross axis coupling in the 
controller as possible, the desired eigenvectors in the longitudinal channel 
may be chosen so that the responses from lateral channels are blocked, i.e., 
select eigenvectors as 


U.I-A) 


C 


LAT 



n 

I 

i=1 


Vi 
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V d are the desired eigenvectors, or are selected using the same criteria as 
in the minimally structure case. Though these two approaches could be made 
to work at each flight condition by iterating on the desired eigenvalues, 
they failed to easily produce a set of acceptable eigenvectors which could 
be used at other flight conditions. Next, partial specification of desired 
eigenvectors was attempted using the dominant mode approach. 

■4 . 2. 1 . 3 - Dominant Mode Eigenassi gnment 

According to ref C 1 3D » complete specification of desired eigenvectors 
are neither necessary nor desirable. Depending on the states that should or 
should not participate in a given mode, appropriate entries in the 
eigenvectors are made ones or zeros, leaving other entries free. With these 
eigenvectors, the desired eigenvalues are moved as far left from the 
imaginary axis as possible with least change in the location of unplaced 
eigenvalues. The desired eigenvectors so obtained appeared to work over 
most flight conditions. Note, however, that extensive iterations on the 
desired eigenvalues may often be required to produce a satisfactory design. 


2. 1.4 Conclusions on Output Eigenstructure Assignment 

Specific maneuver autopilot design results are discussed in Appendix 
A, along with the more complete description and evaluation of output 
eigenstructure assignment for application to flight test trajectory control, 
including difficulties in selecting desirable eigenvalues and eigenvectors. 
This approach demands several iterations to converge to a satisfactory 
design and does not appear to easily give suitable insight for output feed- 
back design of high order multivariable systems which will be used at other 
operating points. If a rational method to generate an achievable set of 
eigenvectors is devised, this technique will be made more attractive. One 
possibility might be to generate gradients of the eigensystem between flight 
conditions and include this information in the single point design techni- 
que. The next two subsections describe two other techniques of output 
feedback gain solution. 
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4.2. 2 Minimum Error Excitation Output Feedback Design 


Following Kosut's [8] notation in describing his development of the 
minimum error excitation output feedback design method, for the linear 
system 


x = Ax + Bu, 
y = Hx, 


one first designs a full state LQ regulator 


if 

2 \ 


(xH T Q Hx 

y 


+ u Ru)dt , 


with the optimal feedback control law u* = F*x. 

Both because algebraic output feedback methods do not guarantee 
stability and because of the presence of slow closed loop full-state 
feedback modes (either neutrally stable unobservable modes or other sluggish 
phugoid like, spiral or integral error modes), it is desirable to design 
with a guaranteed stability margin [14, 15]. Specifying a stability margin 
of a ensures that the real parts of all closed loop roots are less than -a. 
This is equivalent to optimizing the performance index 


j = 1 j e 2at (xH T Q Hx + u T Ru)dt, 

2 o y 
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and can be accomplished by destabilizing the open loop plant by A = A + al; 
solving for the corresponding optimal gains in the standard LQ problem, and 
using them with the original open loop dynamics matrix A. Costs were chosen 
according to Bryson's rule [15], the inverse of the squared deviation 
desired on inputs and outputs, with a scalar factor between Q and R to 
regulate the extent of high gain solution achieved. 

A minimum norm output error feedback law can be determined from 


u = C y, with 

y 

C = H + F* , 

y 


+ 

where H is the pseudoinverse of H. This is merely a least squares projec- 
tion of the full state gains onto the output subspace. 

The sensitivity of the projection described above can be minimized by 
doing a weighted least squares, where the weight is the closed loop state 
covariance excited by a unit error density. One merely solves the Lyapunov 
equation 


(A + BF* ) P + p ( A + BF*) T + 1 = 0, 


to obtain the constrained gain 


F y = F«PH T (HPH T ) -1 H, or 
C - F»PH T (HPH T ) _1 

y 


The destabilized open loop plant A was used instead of A in the Lyapunov 
equation solution to retain as much of the full-state law guaranteed 
stability margin as possible. 
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Maneuver autopilots were designed at all flight conditions within two 
iterations using the above output feedback design procedure. The design 
criterion were 


Re(A^) <-.2 => a 5 second maximum time constant in tracking, and 
£ > .7 => low overshoot. 


With the integral error states and general low frequency oscillatory poles 
pairs due to lateral/longitudinal coupling, the plant is not a simple one to 
control. Since Mach and load factor were chosen as the scheduling 
variables, 15 designs were generated - five designs at three load factors 
shown in Table 4-2. The state and control weighting matrices used in this 
synthesis is listed in Appendix E. The corresponding output feedback gains 
are also given in this appendix. 

Since the zoom-and-push-over and excess thrust windup turn require 
fixed throttle, 15 more designs were generated without this control, also 
shown in Table 4-2. 

Typical initial condition responses in a straight and level flight 
condition at 40K feet and Mach of 1.4 are shown in Figure 4-4. 
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TABLE 4-2. SLOWEST MODE BEHAVIOR AT ALL DESIGN CONDITIONS 

(Subscript 1 corresponds to the design with all con- 
trols, subscript 2 corresponds to the design without 
throttle control) 


h ,ro 

LOAD 

FACTOR 


h = 10000' 
M = 0.8 


h = 20000' 
M = 1.0 


h = 30000' 
M = 1.2 


h = 40000 ' 
M = 1 .4 


h = 50000 ' 
M = 1.8 


1 


X, = -0.277 
Ei = 0.996 

X 2 = -0.28 
E 2 » 1 


A j = -0.289 
Ei = 0.97 

X 2 = -0.257 
E = 1 


X i = -0.3 
Ei = 0.95 

X 2 = -0.234 
E 2 = 1 


Ai = -0.308 

E, = 0.945 

E 2 - -0.2276 
S 2 = 1 


X, = -0.292 
E, = 0.94 

X 2 « -0.219 
E 2 = 1 


2 


X! = -0.28 
Ei = 0.99 

X 2 = -0.288 
E 2 - 0.99 


X , = -0.26 
Ei = 0.94 

X 2 = -0.268 
E 2 - 0.947 


X, = -0.23 
E, = 0.9 

X 2 = -0.254 
E 2 = 0.947 


Xj - -0.2274 
Ei = 0.2 

X 2 * -0.276 
E 2 * 0.958 


X, = -0.225 
X! = -0.85: 

X 2 = -0.242 
E 2 = 0.88 


4 


Xj - -0.267 
Ei = 0.99 

X 2 * -0.28 
E 2 - 0.99 


X, = -0.23 
Ei = 0.92 

X 2 = -0.23 
E 2 = 1 


X, = -0.209 

Ei = 0.866 

X 2 = 0.243 
E 2 = 0.95 


Xj * 0.17 
Ei - 0.98 

X 2 = -0.2405 
C 2 * 1 


X, * -0.116 
Ei = 0.95 

X 2 - -0.209 
E 2 * 0.895 
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Figure 4-4. Initial Condition Response for the Output Feedback Minimum 
Error Excitation Perturbation Controller (Reference 
Condition H = 40k', Mach = 1.4, Straight and Level Flight) 


Section 5 gives the tracking behavior of the output feedback minimum 
norm controllers in the eight maneuvers of interest. 


4.3 NONLINEAR FLIGHT TEST TRAJECTORY CONTROLLERS 

Research on nonlinear flight test trajectory control (FTTC) design was 
conducted with the goal of completely eliminating the need for gain 
scheduling. A brief literature review is given below, outlining how recent 
theoretical work can be applied to the FTTC problem before a demonstration 
problem is given to illustrate the approach. 
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This research deals with the synthesis of nonlinear flight test 
trajectory controllers using the recent results in prelinearizing transforms 
due to Meyer [ 1 7 — 333 and singular perturbation theory [3^ _ 36], The use of 
singular perturbation theory in this problem simplifies the command 
generation scheme in addition to providing a consistant approach for 
eliminating ignorable state variables. The prelinearizing transformations 
are more transparent in this formulation. The slow-fast computations are 
clearly separated and can be carried out at different rates on the flight 
control computer. It is interesting to note that in Ref. 21, even though 
the controller development did not make use of singular perturbation theory, 
the time-scale separation formed a basis for implementation on the flight 
control computer. A schematic block diagram of the slow-fast flight test 
trajectory controller is given in Figure H-5. 

Flight test controller synthesis will be developed for the F-15 
fighter aircraft in the next contract phase. For the F-15 and most fighter 
aircraft, it can be assumed that the aircraft under consideration has the 
four usual controls: throttle, aileron, rudder and elevator. The objective 

of the flight test controller is to track the given commands in airspeed, 
angle of attack, angle of sideslip and altitude in presence of disturbances 
and modeling imperfections. It is clear that the commanded trajectory has 
to be executable by the aircraft under consideration. Note that the flight 
test control problem discussed here is distinct from those described by 
Meyer, et . al . [20-22] since, in their work, the trajectory to be followed 
consisted of the three position components specified as functions of time. 

Modeling and time-scale separation can be exploited to a considerable 
degree in this approach. The mechanization details of prelinearization and 
the slow-fast controller synthesis for general flight test maneuvers will be 
given in the next project phase based on the problem formulation developed 
here. 



Figure H~5. Singular Perturbation Nonlinear Flight Test Trajectory 
Controller . 

Slow states: V, a, 8, h 

Fast states: p, q, r, 0, <|> 

Very fast states: CAS and Actuator states 


An illustrative example of this nonlinear controller approach is given 


in Appendix B. 







SECTION 5 

SIMULATION AND EVALUATION 


The performance of perturbation controllers for initial condition 
errors was verified during the design phase. Some techniques for evaluating 
their tracking performance will be discussed in this section. Though it is 
clear that the actual performance of these controllers can only be assessed 
with the full nonlinear aircraft model, it is desirable to introduce an 
intermediate validation phase to ensure a smooth transition from 
perturbation controller design to full nonlinear simulation. It should be 
emphasized at the outset that the simulations discussed here are approximate 
and consequently the tracking performance will differ in the nonlinear 
simulation validation in the next study phase. 


5.1 MANEUVER SIMULATION 


Originally, it was decided to carry out a linear time varying 
simulation of the systems using linearized aircraft models along the desired 
flight-test trajectory as discussed in Appendix D. For a two state variable 
model analyzed in Appendix B, the computing time was small. However, the 
aircraft and CAS system has 31 states and continuous interpolation was found 
to be extremely time consuming. . In view of the excessive computational 
effort required, and the limited value of the information obtained, it was 
then decided to switch models and controller gains along a desired flight 
test trajectory. This approach introduced artificial gain switching 
transients and also led to misleading conclusions. Hence whenever feasible, 
the simulations discussed here used one interpolated model, and gains based 
on the flight condition halfway through the maneuver. With this approach, 
the modeling inaccuracies will be almost equally distributed throughout the 
maneuver . 
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It should be emphasized that the simulations with linearized models 
will test only the feedback controller portion of the maneuver autopilot. 
The open loop control histories generated from maneuver modeling can be 
tested only during the full nonlinear simulation of aircraft and CAS. 


5.2 MANEUVER SIMULATION MECHANIZATION 


The linear perturbation equations used for design have been discussed 
in Section H. In this section, these will be modified to generate linear 
simulations. The linearized aircraft with CAS is of the form 


Sx ■ F6x + G6u 
6y * H6x 


The output feedback perturbation controller is of the form 


6u = C y 6y. 


Expanding the perturbation equations back out gives 


x - x * F(x-x) - GC y H(x-x) 

X = (F - GC y H)x - [(F-GC y H)x - x] , 
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from which the maneuver simulations can be mechanized. Note that x can be 
viewed as a disturbance (and neglected) or as a part of the external 
reference command. For the high speed maneuvers simulated here x clearly 
cannot be ignored, but were computed numerically with a forward differencing 
of x. 


5.3 MANEUVER SIMULATION RESULTS 

As discussed earlier, x was computed numerically from x and used as 
part of the reference command. Since x has corners, there are jumps, or 
spurious step inputs due to the discontinuities in x. Using a smooth x from 
quadratic or cubic spline fits to the trim points will remedy this 
problem; however, the effects of the discontinuities can clearly be seen in 
the tracking trajectories in this subsection. 

In the following, typical simulation results for each flight test 
trajectory will be presented. The flight conditions in these simulations 
are chosen so that the maneuver autopilot is exercised over nearly the 
entire aircraft envelope. Except where indicated, in the plots that follow, 
the dotted lines denote the commanded variables generated from the maneuver 
modeling program while the solid line represents the trajectory evolution 
from the linear simulation. 


5.3.1 Transient Trajectory 

In this simulation, a transient trajectory is setup to transfer the 
aircraft from straight and level flight conditions at 20000’ altitude and 
Mach 0.8 to straight and level flight conditions at 30000' altitude and Mach 
1.2 in 60 seconds. The simulation results are presented in Figs. 5~1 
through 5-3. The altitude tracking is very good. However, the throttle 
history in Fig. 5-3 indicates that during the first five seconds, the 
controller commanded a negative throttle. This is caused by the nonminimum 
phase behavior of the aircraft, clearly discernable in the Mach number and 
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angle of attack histories. In the actual implementation, the commanded 
altitude can be modified to include an initial descend leg in order to 
compensate for the nonminimum phase behavior. Alternately, the commanded 
Mach number can be modified to have an initial decreasing segment. In any 
case, mere command modification would fix the initial negative throttle 
difficulty. 
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TIME, SECONDS 



Figure 5 - 1. Altitude and Mach Number Evolution Along the Transient 
Trajectory 
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Fig U*"' a 5~2. Angle of Attack and Flight Patn Angle Evolution Along the 
Transient Trajectory 
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Figure 5-3. Throttle and Elevator Deflection Along the Transient 
Trajectory . 

5.3.2 Level Acceleration 

A level acceleration flight test trajectory at 30000 1 altitude is 
considered here. The aircraft is required to accelerate from Mach 0.9 to 
Mach 1.2 in 60 seconds. The tolerence on altitude is +_ 50' while Mach 
number error should be within ^+0.01 . The simulation results for this 
maneuver are presented in Figs. 5"^ through 5-6. The Controller was able to 
maintain the altitude within +0.1 feet while tracking the Mach number within 
the given specifications. The throttle and elevator deflections given in 
Fig. 5-6 are well within the saturation levels. From the maximum throttle 
requirment in this maneuver, it appears that the maneuver time could be 
decreased by about 20 seconds. 
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Figure 5~6. Throttle and Elevator Deflection Along the Level 
Acceleration Flight Test Trajectory* 


5.3.3 Pushover, Pullup 

This maneuver was initiated at 30000* altitude and 0.8 Mach straight 
and level flight condition. The objective is to track a piewise linear 
angle of attack history given in Fig. 5-8 while maintaining the Mach number 
constant at the initial value. The results of the maneuver simulation are 
given in Figs. 5-7 through 5~9. From Fig. 5~7, it can be seen that the Mach 
number error is within 0.001 of the commanded value. The angle of attack 
tracking error is less than 0.2° throughout the maneuver. The throttle 
history given in Fig. 5~9 shows a small negative region at the imitation of 
pullup at 30 seconds and is primarily due to the corner in the angle of 
attack command. Smoothing this corner in the final mechanization would 
eliminate this negative throttle region. 
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Figure 5“9. Throttle and Elevator Deflection along the Pushover/Pullup 
Flight Test Trajectory. 


5.3-^ Zoom and Pushover 

As noted in the maneuver modeling, this flight test trajectory is 
executed in three phases. In the first phase the aircraft is transferred 
from straight and level flight conditions to the beginning of the zoom and 
pushover parabolic trajectory. The second phase consists of the zoom and 
pushover trajectory while the third phase restores the aircraft to the 
original straight and level flight condition. During the first and third 
phases all the controls are active, but the throttle is fixed during the 
zoom and pushover trajectory. 
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In the present case, the aircraft is first trimmed to fly straight and 
level at 30000' altitude and Mach 0.4. An initial transient trajectory is 
then executed with all controls active until the point marked A in Fig. 5 - 
10. At this point, the throttle is fixed and the aircraft executes the 
parabolic zoom and pushover trajectory until the point B in this figure. 

The throttle is released at point B and the aircraft performs another 
transient maneuver to restore it to the original conditions. 

Since the controller performance along the transient trajectory has 
already been investigated in Section 5.3.1, the tracking performance along 
the zoom and pushover trajectory will only be demonstrated here. The 
aircraft begins the zoom and pushover maneuver at about 28000' and Mach 0.45 
and completes the maneuver at about the same conditions. The controller 
performance is illustrated in Figs. 5 _ 10 through 5-12. The conditions at 
the apex of the parabola is of particular interest in this maneuver. From 
Fig 5-11, it can be observed that the angle of attack at the apex is within 
0.005 radian of the required value. The altitude error is within 50' at the 
apex and the Mach number is within 0.05 of the required value. Initial 
transients in altitude and Mach number can be seen at point A in Fig. 5-10. 
These are essentially due to the availability of just one control variable, 
the elevator, to track three state variables: altitude, Mach number and 

angle of attack. Thus, an initial condition error on altitude would 
propagate to Mach number and angle of attack channels and vice versa. 
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Figure 5-12. Elevator Deflection along the Zoom and Pushover 
Flight Test Trajectory. 

5.3.5 Excess Thrust Windup turn 

The results for an excess thrust windup turn trajectory at 40000' 
altitude and Mach 1.4 are given in Figs. 5-13 through 5-19. The altitude 
and Mach number were required to be constant throughout the windup turn 
trajectory, while tracking an angle of attack command as shown in Fig. 5-14. 
The aircraft roll altitude in this maneuver is close to 70° and results in a 
highly coupled model to be controlled by the maneuver autopilot. It can be 
observed that the maneuver autopilot maintained the altitude within +5' and 
the Mach number error is within +0.002. Except at the beginning and the end 
of the maneuver, the angle of attack error is within 0.005 radian of the 
commanded value. For about 10 seconds during the beginning and end of the 
maneuver, the rudder deflection is close to 12° while in the high roll 
attitude region, an elevator deflection of nearly 17° was demanded. This 
indicates that at the present flight conditions, a less stringent maneuver 
should be attempted. 
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Figure 5-13. Altitude and Mach Number Eve- ut ion along the Excess Thrust 
Windup Turn Flight Test Trajectory 
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Figure 5-lH Angle of Attack and Flight Faun Ang-e Evolution along the 

" ' Excess Thrust Windup Turn Flight Test Trajectory 
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Figure 5—15. Roll Attitude and Angle of Side Slip Along the Excess Thrust 
Windup Turn Flight Test Trajectory 
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Figure 5-1 6. Roll and Pitch Body Rate Evolution A-tong the Excess Thrust 
Windup Turn Flight Test Trajectory 
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Figure 5 — 17 . Yaw Body Rate Evolution along the Excess Thrust Windup Turn 
Flight Test Trajectory 
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Figure 5-18. Throttle and Elevator Deflection along the Excess Thrust 
Windup Turn Flight Test Trajectory 
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Figure 5*19. Differential Tail and Rudder Deflection along the 
Excess Thrust Windup Turn Flight Test Trajectory. 

5.3.6 Constant thrust windup turn 

A descending constant thrust windup trajectory will be illustrated in 
the following. The aircraft starting at straight and level flight condition 
enters a level turn with linearly increasing angle of attack, upto 30 
seconds in Fig. 5*21. At this point, the throttle is fixed and the constant 
thrust windup trajectory begins. At the end of the maneuver , the Angle of 
attack is gradually decreased to the straight and level trim values. During 
the constant throttle windup turn trajectory, the Mach number is to remain 
constant. The simulation results for this maneuver are given in Figs. 5*20 
through 5*26. From Fig. 5*20, it can be seen that the Mach number was 
maintained within +0.0075 while the angle of attack tracking error was 
within 0.01 radians. The control surface deflections were within the 
saturation limits except at the point where the constant throttle windup 
trajectory began. This is due to the corner present in the altitude and 
angle of attack command histories. By smoothing these corners in the 
commands, the control surface limit violation can be avoided. 
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Figure 5“20. Altitude and Mach Number Evolution Along the Constant 
Throttle Windup Turn Flight Test Trajectory. 
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Figure 5~21 . 


Angle of Attack and Flight Path Angle Evolution Along the 
Constant Throttle Windup Turn Flight Test Trajectory. 
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Figure 5-22. Roll Altitude and Angle of Side Slip Evolution Along the 
Constant Throttle Windup Turn Flight Test Trajectory. 
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Throttle Windup Turn Flight Test Trajectory. 
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Yaw Body Rate Evolution along the Constant Throttle 
Windup Turn Flight Test Trajectory. 
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Figure 5-25. Throttle and Elevator Deflection along the Constant Throttle 
Windup Turn Flight Test Trajectory. 
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Figure 5-26. Differential Tail and Rudder Deflection along the Constant 
Throttle Windup Turn Flight Test Trajectory. 


5. 3.7 Constant Dynamic Pressure - Constant Load Factor Maneuver 

As noted in the maneuver modeling section, this maneuver can be 
ascending or descending based on the required Mach rate. The controller 
performance along an ascending constant dynamic pressure-constant load 
factor trajectory is given in Figs. 5~27 through 5~34. This trajectory 
consists of three phases. The first phase begins with initial conditions 
chosen to obtain the desired dynamic pressure. Next the aircraft is placed 
in a turn to generate the required load factor. In the present case, a load 
factor of 4 was employed. Next, the desired Mach rate is initiated and the 
constant dynamic pressure-constant load factor trajectory is executed. In 
the present case, in order to achieve a Mach rate of 0.0067/second, the 
aircraft had to climb from 35000' altitude to 43000' altitude in 30 seconds. 
These histories are given in Fig. 5-27. The control surface deflections are 
well within the saturation limits except at the points where the altitude 
and Mach number commands contains corners. The dynamic pressure history 
given in Fig. 5~34 was maintained within 4$ of the required value. 
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e 5-31. Yaw Body Rate along the Constant Dynamic Pressure Constant 
Load Factor Flight Test Trajectory. 
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Figure 5 - 32 . Throttle and Elevator Deflection Along the Constant 
Dynamic Pressure Constant Load Factor Flight Test 
Trajectory. 
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Dynamic Pressure Constant Load Factor Flight Test Trajectory 
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Figure 5~3^ . Dynamic Pressure along the Constant Dynamic Pressure Constant 
Load Factor Flight Test Trajectory 
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5.3.8 Constant Reynold's Number - Constant Load Factor Maneuver 


This maneuver sequence is identical to the constant dynamic pressure- 


constant load factor trajectory. The controller performance for a negative 
Mach rate - Constant Reynolds' number-constant load factor trajectory is 
given in Figs. 5~35 through 5-42. Note that the Reynold's number given in 
Fig. 5-42 should be multiplied by the characteristi c diameter to obtain the 
actual Reynold's number. From the performance results given in Fig. 5-40 
and 5-41, it can be seen that the control surface deflections momentarily 
violate the saturation limits at the points corresponding to the corners in 
the commanded values. Throughout the trajectory, Reynold's number is within 
1.5$ of the required value. 
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Figure 5 - 35 . Altitude and Mach Number Evolution Along the Constant 
Reynold's Number Constant Load Factor Flight Test 
Trajectory. 


82 




BETA , RAB PHI , RAO * GAMMA , RAD ALFA , rad 


Figure 



0 10 20 30 40 50 60 70 BO 90 100 



0 10 20 30 40 50 60 70 BO BO 100 

7ItfP , SECONDS 


'“36. Angle of Attack and Flight Path Angle Evolution along the 
Constant Reynold's Number Constant Load Factor Flight Test 
Trajectory 
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-37. Roll Attitude and Angle of Side Slip Evolution Along the 

Constant Reynold's Number Constant Load Factor Flight Test 
Trajectory . 
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Figure 5-38. Roll and Pitch Body Rates along the Constant Reynold's Number 
Constant Load Factor Flight Test Trajectory. 



Figure 5~39. Yaw Body Rate along the Constant Reynold's Number Constant 
Load Factor Flight Test Trajectory. 
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SECTION 6 

SUMMARY AND FUTURE WORK 


This report summarizes the analysis done in the first phase of a study 
for developing and validating a flight test trajectory controller for eight 
maneuvers. System modeling, control design results, control design 
technique evaluation, validation and software deliverable results and 
conclusions are summarized here. This final section describes future work 
which will build upon this study in the further validation and flight 
testing of these flight test trajectory controllers. 


6.1 SYSTEM MODELING 


With newly developed linearization tools by NASA Ames Dryden Flight 
Research Facility, the linear airframe models are straight forwared to 
obtain. Command augmentation system CAS models are important and although 
they complicate the maneuver autopilot designs, they are essential since the 
maneuver autopilot must work through the CAS in the manned simulation or in 
the F-1 5 itself. 

The output analytical constrained relations developed in section 3 
which discusses maneuver modeling were found to be fundamental in achieving 
a consistent maneuver autopilot. The nonlinear aerodynamics can be 
effectively reduced to a reference command table by appropriate 
linearizations throughout the flight envelope. A small number of linear 
perturbation models can be used to decompose the eight desired maneuvers 
into a few sets of linear perturbation equations. 


6.2 CONTROL DESIGN RESULTS 


Output feedback controllers with appropriate integral aerostates were 
found to be the simplest form of feedback controllers. Integral aerostates 
must be chosen carefully to avoid control lability problems which are 
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problematic in the design stage. The use of a guaranteed stability margin 
in the control design technique was found to be a powerful method when 
working with the high order augmented and coupled models. It was found that 
a five second response time with a minimal amount of overshoot was 
straightforward to achieve at all design points. 


6.3 CONTROL DESIGN TECHNIQUE EVALUATIONS 


The eigenstructure assignment technique was found to work well on a 
bare airframe but with a complex augmented model was not useful as a design 
technique because there is no procedure to perform a converging iterative 
design approach. 

The minimum error excitation output feedback method with a preliminary 
guaranteed stability margin full state design worked well at all flight 
conditions. 


6.4 MANEUVER AUTOPILOT VALIDATION 


The maneuver autopilots were validated in a linear simulation. The 
tracking response for reasonably high performance maneuvers were found to be 
quite acceptable, and except for exceeding the control authority at isolated 
points, is ready for further testing in a full nonlinear simulation. 


6.5 SOFTWARE DELIVERABLES 


The tools used in model development control design evaluation and 
design throughout the envelope as well as validation were made available to 
NASA Ames Dryden Flight Research Facility. These included a time varying 
simulation mechanized in ISI's MATRIX^, SYSTEM_BUILD (see Appendix A). In 
addition a stand alone program which performs a three dimensional 
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interpolation in state variables and converts this through the interpolation 
into a ID table as a function of time for a specific maneuvers was developed 
and also delivered to NASA. The 1-D interpolation in time can then be used 
in SYSTEM_BUILD to mechanize a time varying simulation. Documented command 
files were delivered for the model generation, control law design, and 
construction of linear simulation using a single model throughout the flight 
maneuver. Command files for building, loading the data and executing a time 
varying simulation model were also provided. 


6.6 FUTURE WORK 


The next phase of this study will involve validation of flight test 
trajectory control laws in a batch nonlinear simulation. The gain scheduled 
linear perturbation controllers developed in this study Are currently ready 
for validation in such a simulation. The fineness of the discretization 
both of the reference commands and linear perturbation models may need to be 
revised in the next phase. An open-loop simulation of the reference control 
values will give an initial check of the accuracy of our nonlinear tabular 
model. The nonlinear simulation provided by NASA Ames Dryden Flight 
Research Facility will include both the CAS model as well as the airframe 
dynamics and Sensor models. The structure of the nonlinear control law, 
while only outlined here, was developed to a point where with some symbolic 
or algebraic manipulation by hand, it can be mechanized on the nonlinear 
simulation. The advantage of this control law is that a single set of gains 
will work for all maneuvers throughout the envelope. Such a transformed 
linear system can be easily controlled and the resulting gains back 
transformed with the nonlinear equations to give a nonlinear control law. 
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APPENDIX A 


FLIGHT TEST TRAJECTORY CONTROLLER SYNTHESIS 
WITH CONSTRAINED EIGENSTRUCTURE ASSIGNMENT 


A-l 


INTRODUCTION: 


This Appendix gives the details of research conducted on the 
eigenstructure assignment. The material in the Appendix A-I was presented 
as a paper in the 1985 American Control Conference at Boston. Appendix A-II 
and A-III give the aircraft - CAS models and the constrained eigenstructure 
designs at two flight conditions. Appendix A-IV gives the desired 
eigenvalues and eigenvectors used in the synthesis. 

The aircraft - CAS model used here is of the form 


x - Fx + Gu 
y = Hx 

where: 

the state variables x are the perturbed values of 

V, total airspeed 

а. Angle of attack 
q, pitch body rate 

б, pitch attitude 

B , angle of sideslip 
p, roll body rate 
Y , yaw body rate 
4>, roll attitude 
h, altitude 

engine actuator state and 21 CAS states 

U: the control vector consists of perturbed values of throttle, 
elevator, differential tail and rudder. 


y: the measurement vector and consists of the perturbed values 

of 

p: roll angular acceleration 

A : normal acceleration 

n 

q: pitch rate 

q: pitch angular acceleration 

p: roll rate 

A y-body acceleration, not at the vehicle center of 

y»i 

gravity 

r: yaw angular acceleration 

r: yaw rate 

h: altitude 

M: mach number 

a : angle of attack 

n: load factor 

<p: roll attitude 

6: pitch attitude 

h: altitude rate 

U3: x-body axis velocity 

A : x-body axis acceleration at vehicle C.G. 

nx J 


CONCLUSIONS: 


The constrained ei genstructure assignment design procedure in its 
present form demands several iterations to converge to a satisfactory design 
and does not appear to easily yield suitable insight for output feedback 
design of high order multivariable systems. If a rational method to 
generate an achievable set of desired eigenvectors is devised, this 
technique will be made more attractive. 
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ABSTRACT 


Flight test trajectory controller synthesis using constrained eigenvalue/ 
eigenvector assignment procedure is presented. Associated modeling problems 
and the difficulties encountered in employing this synthesis technique are 
highlighted. 
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1. INTRODUCTION: 


Flight test trajectory control is a technique designed to aid in the 
collection of large quantities of high quality data. This technique has 
provided the means for flying maneuvers consistantly , precisely, and repeat- 
ably from flight to flight. Two versions of these controllers have been 
used: a closed-loop automatic system and an open-loop system providing 

manual piloting information. A closed-loop system used to collect per- 
formance, pressures, and loads data from the highly maneuverable aircraft 
technology (HiMAT) vehicle is described in [1], The application of the 
open-loop system on the NASA F— 111 transonic aircraft technology (TACT), 

F-15 airframe/propulsion system interaction studies, and F-15 shuttle tiles 
test programs are given in [2], 

Originally, the open-loop flight-test- trajectory guidance algorithms 
were developed on-line, in a piloted simulation using cut-and-try techniques 
that was not only man power intensive, but often produced less than desir- 
able controllers. A closed-loop system designed using one-loop-at-a-time 
classical design approach is documented in [3]. Full-state feedback 
approach for closed-loop system design using Linear quadratic synthesis is 
described in [A]. Both these approaches have limitations in terms of design 
methodology and controller complexity. 

The research currently underway includes an exploration of various 
multivariable synthesis techniques for this problem. The first approach 
considered is that of constrained eigen value/eigen vector assignment [5], 

A primary goal is to develop controllers based on output feedback so as to 
decrease controller complexity and to enhance robustness. The objective 
of this paper is to point out strengths and weakness of this technique when 
used in a, realistic, relatively a large problem such as required for flight 
test trajectory controllers. 
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2. THE FLIGHT TEST TRAJECTORY CONTROL PROBLEM : 

Flight test trajectories are flown to evaluate an aircraft within its 
known operational envelope and to explore the boundaries of its capabili- 
ties. This makes the flight test trajectory control a very demanding task. 
Control systems designed for this purpose must not only operate satisfac- 
torily in terms of keeping the flight test variables within acceptable 
tolerance but should also be reasonably insensitive to model parameter 
variations . 

An approach to closed-loop flight test trajectory controller synthesis 
consists of linearizing the aircraft model at several flight conditions 
about the flight test trajectory and designing multivariable controllers, 
which in some sense minimizes the deviations from the reference path. Time 
varying nature of the linearized aircraft model along the reference path 
brings about the need for scheduling gains as a function of time or as func- 
tions of some important flight variables such as dynamic pressure, Mach 
number etc. The gain scheduling aspects will not be pursued any further in 
this paper, and in all that follows, discussions will center around a 
linear-time-invariant aircraft model. Further, though the flight test tra- 
jectory controller is discrete, in this preliminary stage it will be assumed 
that the sampling rate is sufficiently high, permitting the application of 
continuous control design techniques. 

The aircraft under consideration is a high performance fighter with 
command augmentation system (CAS) engaged in all the three axes. The CAS is 
a highly nonlinear system with saturations, multiplicative nonlinearities 
and gain schedules. At a particular flight condition, this system can be 
approximated by a linear system with "equivalent" gains derived from non- 
linear simulations. The aircraft model at the same flight condition is 
obtained from a generic aircraft linearization code developed at NASA-Dryden 
Flight Research Facility [6]. The state variables used are total speed, 
angle of attack and angle of sideslip, pitch rate, yaw rate, roll rate, 
pitch attitude, roll altitude and altitude. Throttle, rudder, elevator and 
differential tail constitute the control variables. The engine dynamics are 
modeled as a first order lag. Thus, the aircraft together with CAS at a 
particular flight condition is a coupled 31 -st order system with four con- 
trols, the aircraft model having 10 states. A block diagram of the system 
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I with flight test trajectory controller (called the Maneuver Autopilot here) 

included in the loop is shown in Figure 1. Note that MAP controls the air- 
craft through CAS. Due to this, even if the aircraft model is decoupled in 

i 

longitudinal and lateral axes at a particular flight condition, the CAS in- 
troduces a strong coupling into the system. In concise terms, the chief 
objective of the present study is to synthesize the maneuver autopilot so as 
to effect satisfactory transient response. Six maneuvers have been 
analyzed, which are sketched below. 

2.1 Level Acceleration/Deceleration 

This is a wings-level, constant altitude maneuver with Mach number 
constant or changing at a specified rate. 

2.2 Pushover, Pullup 

This is a wings-level, constant Mach number maneuver in which angle of 
attack is varied a specified increment about the trim value at some speci- 
fied rate. 

2.3 Zoom and Pushover 

The zoom and pushover is a wings-level, thrust stabilized less than lg 
maneuver. The flight trajectory is a parabolic path with the target 
Mach/altitude/angle of attack point at the apex. 

1 2.4 Excess Thrust Windup Turn 

i 

j 

This is a maneuver with angle of attack linearly increasing from the 
wings-level trim condition to some specified final value at a specified 
rate. The maneuver is performed at constant altitude and constant Mach 
number. 
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2.5 Constant Throttle Windup Turn 


This is a maneuver with angle of attack increasing linearly at a 
specified rate from trim to some specified final value. The maneuver is 
performed at a predetermined, constant thrust level. Mach number is main- 
tained by trading potential for kinetic energy via an appropriate altitude 
rate. 

2.6 Constant Reynolds Number and Constant Load Factor Trajectory 

This maneuver is initiated at a predetermined load factor, Mach num- 
ber, and dynamic pressure. Thus, the initiation of this maneuver is not 
necessarily the wings-level condition. This maneuver can be either an 
ascending or descending maneuver at a specified Mach number rate. Reynolds 
number and load factor are held constant throughout the maneuver. Altitude 
is gained or lost to maintain Reynolds number with changing Mach number. 

The simplest of these is the level acceleration/deceleration trajec- 
tory and the MAP design for this maneuver will be used as the illustrative 
example in this paper. 


3. OUTPUT FEEDBACK DESIGN: 


As indicated earlier, output feedback is attractive because of simpli- 
city. Further, gain scheduling will be essential in this situation and in 
order to minimize the amount of stored data, it is desirable to have a 
capability to impose control structural constraints. 

Currently, several approaches are available for output feedback 
design, see Refs [5, 7~ 1 0] for example. It is not the purpose of this paper 
to compare and contrast these, but to evaluate a specific technique from an 
application point of view. Constrained eigen value/eigen vector assignment 
technique of [5] was used in this particular application with the hope that 
it would permit the generation of designs based on the classical notions of 
poles, zeros and their relative location in the complex plane. 

-Additionally, the selection of this technique was motivated by the example 
in [53, viz, the lateral axis stability Augmentation System for the L1011 
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airplane. Salient features of this approach are outlined in the following 
for clarity of the subsequent discussions. 

3.1 Constrained Eigen Value Vector Assignment [5]: 

Consider the linear time invariant system 

x = Ax + BU 

y « Cx 


with X e R n , U e R m , y e R r and A, B, C are real constant matrices of 
compatible dimensions. 

It is desired to design a feedback controller of the form 


U - Fy 


with a structural constraint that some specified elements of F satisfy 


f 


U 


0 


It is assumed that the system can be stabilized with given outputs y, 
and that any dynamic compensators required have been appended to the 
original system. 

The design problem is: given a self conjugate set of scalars 

{A^} i « 1, 2 . ..r, and a corresponding self conjugate set of n vectors 
{vp i « 1, 2 ... r, and a given controller structural constraint, find a 
real m x r matrix F such that r of the eigen values of A + BFC are 
"close” to the set { A^} and the corresponding eigen vectors of A + BFC are 
"close" to the respective member of {vp. 
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Note that if there are no constraints on the feedback matrix, it is 
feasible to place r desired eigen values exactly. 

However, one is not at liberty to place v , v ... v by an arbitrary 

d 1 t ^ 

set of desired closed loop eigen vectors {v..} i = 1, 2 ... r since this set 
might not belong to the set of assignable closed loop eigenvectors. It then 
becomes necessary to find some approximation to the set {v^} which is 
assignable, yet "close" in some sense to the original {v^. It will be seen 
subsequently that this is the main difficulty in the application of this 
technique to a complex problem such as flight test trajectory controller 
design . 

Complete details of this approach can be found in [5]. It is perhaps 
of interest to note that in addition to the capability to handle controller 
structure constraints, this technique can be modified to accept partial 
specification of eigen vectors. This feature can be valuable in high order 
systems. In summary, to employ this synthesis approach, the following are 
required. 

(i) Based on the practical aspects of the problem, choose a minimal 
set of measurements which will permit the designer to achieve the 
desired performance. Introduce dynamic compensators such as 
integral feedbacks, lead-lag networks, etc., based on experience. 

(ii) Choose a set of desired eigenvalues and eigenvectors equal to the 
number of outputs. 

While the selection of desired eigenvalues is often apparent in a 
given problem, the desired eigenvectors are difficult to select. A piece of 
information which might serve as a guide in this selection is the fact the 
closed loop eigenvectors v corresponding to the closed loop eigenvalues A. 

1 -i 1 

must lie in the subspace spanned by the columns of the matrix (A^I-A)B . 
Further, the matrix C V, V = {v^ v^ ... v r } should be invertible. If the 
complete specification of desired eigenvectors are not feasible, a rough 
rule is to pick the entries in these eigenvectors as zeros or ones based on 
whether a particular measurement needs to contribute to a particular mode or 
not. In any case, constructing a set of desired eigenvectors that are 
assignable constitutes the most difficult part of this design technique. 
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1*. FLIGHT TEST TRAJECTORY CONTROLLER SYNTHESIS: 


As an illustrative example, constrained ei genstructure assignment 
technique is next used in the design of level acceleration/deceleration 
flight test trajectory. Alternate design approaches have been used for this 
maneuver and hence comparative evaluation was feasible. The level 
acceleration/ deceleration maneuver requires the roll attitude to be main- 
tained zero and the altitude to be held constant. The mach number should be 
maintained constant or should change at a specified rate. To ensure zero 
steady state errors, integral feedbacks are first introduced in the altitude 
and mach number channels. Since the mach number command can be a ramp, an 
additional integrator is required in this channel to decrease tracking 
errors. This, however, was not done because the tracking errors with a 
single integral feedback has been found to be within acceptable values. 

The desired eigenvalues to be used in the design are clear at the out- 
set, based on the four modes for aircraft model, viz, the phugoid, short 
period, dutch roll and roll convergence. The controller structure con- 
straint in this flight maneuver is that the errors in lateral channel will 
be corrected using rudder and differential tail, while the errors in longi- 
tudinal channel will be corrected using elevator and throttle. In view of 
the time varying nature of the model, one would like to pick a set of eigen- 
values and eigenvectors at a particular flight condition through extensive 
design iterations and then attempt to use these at other flight conditions. 

The choice of desired eigenvectors is not clear at this point. Three 
approaches were tried with varying degrees of success. These are sketched 
in the following. 

ij.1 Minimally Restructered Eigenassignment 

Since it is known that the closed loop eigenvectors lie in a subspace 
spanned by the columns of (AjI-A) B 1 , i - 1, ... n, linear combinations of 
these vectors were used as desired eigen vectors. The weights to be used in 
generating these linear combinations were constructed from the additional 
information that unassigned eigenvectors should be close to their open loop 
values in a least square s6nse. 
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U.2 Decoupling Eigenajssignment 


A variation to the above was attempted next. Since we are interested in 
having the least cross axis coupling in the controller as possible, the 
desired eigenvectors in the longitudinal channel may be chosen to that the 
responses from lateral channels are blocked, i.e. select eigenvectors as 


(X.I-A) J 

B 


" V. ~ 
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0 
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_ C LAT | 
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w . 
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n 

I a.v. 
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are the desired eigenvectors, or are selected using the same criteria as 
4.1. Though these two approaches could be made to work at each flight con- 
dition by iterating on the desired eigen values, they failed to easily pro- 
duce a set of acceptable eigen vectors which could be used at other flight 
conditions. Next, partial specification of desired eigenvectors was 
attempted as follows: 

• 

4.3 Dominant Mode Eigenassignment 

According to ref [5], complete specification of desired eigenvectors are 
neither necessary nor desirable. Depending on the states should or should 
not participate in a given mode, appropriate entries in the eigenvectors are 
made ones or zeros, leaving other entries free. With these eigenvectors, 
the desired eigenvalues are moved as far left from the imaginary axis as 
possible with least change in the location of unplaced eigenvalues. The 
desired eigenvectors so obtained appeared to work over most flight condi- 
tions. Note, however, that extensive iterations on the desired eigenvalues 
may often be required to produce a satisfactory design. 
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The open-loop eigenvalues for the aircraft-CAS linearized model at 
Mach 1.2 and 10000 ' altitude is given in table 1. The system has a pole on 
the right half of s plane. Partial specification of the desired eigen- 
vectors are constructed next, based on the states that should or should not 
participate in a given measurement. The desired eigenvectors used in the 
present example are given in table 2 . It can be seen from this table that 
no restriction has been placed on states associated with CAS. A manual 
iteration is now undertaken to determine the desired eigenvalues. Since the 
designer is interested in producing a stable system with as high a speed of 
response as possible, the eigenvalues to be moved are the ones closest to 
the imaginary axis. Hence, these are moved as far to the left of the 
imaginary axis as possible with least change in the location of other eigen- 
values. Table 3 shows a set of desired eigenvectors obtained from this 
exercise. The output feedback gains obtained from the constrained 
eigenvalue/eigenvector design technique is given in table iJ. Table 5 gives 
the closed loop eigenvalues. Comparing this with table 3 shows that the 
achieved eigen values are close to the desired ones. 

In Figures 2 and 3» the time response of the system for a ramp mach 
number command are shown. A question that occurs naturally at this point is 
whether the design can be improved by further adjustment of desired eigen 
values and eigen vectors. Examination of the constrained eigen value/eigen 
vector approach yields no answer to this question. 

5 . CONCLUSIONS 


Design of a maneuver autopilot for flight test trajectory control 
using constrained eigenvalue/eigenvector assignment was discussed. Diffi- 
culties encountered in the generation of desired eigenvalues and eigen- 
vectors were outlined. This approach demands several iterations to converge 
to a satisfactory design and does not appear to easily give suitable insight 
for output feedback design of high order multivariable systems which will be 
used at other operating points. If a rational method to generate an 
achievable set of eigenvectors is devised, this technique will be made more 
attractive. One possibility might be to generate gradients of the eigen- 
system between flight conditions and include this information in the single 
point design technique. 
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TABLE 1. OPEN-LOOP EIGEN VALUES FOR AIRCRAFT-CAS LINEARIZED MODEL AT MACH 
1 .2 AND 1000’s ALTITUDE 


-5.94560+01 +2.0293D+01i 
-5.94560+01 -2.02930+0ii 
-1.36740+02 -3. 262SD— 15i 
-4.88790+01 +3.04340-17i 
-8.77580+01 +5.82180+0H 
-8.77580+01 -5.82180+0U 
-1.01800+02 +5.17230-15i 
-1.00000+02 +4.21440-15i 
-3.81980*01 +3.54690-16i 
-2.32860+01 +8.24360-17i 
-7.41170+00 +3.35740+01i 
-7.41170+00 -3. 35740+01 i 
-7.97250+00 +8.2€550+00i 
-7.97250+00 -8.26550+001 


-1.21760+01 +1.01920+01i 
-1.21760+01 -1.01920+01i 
-2.94380+00 +5.31130+001 
-2.94380+00 -5.31130+00i 
-2.31860-01 +1. 85040-4)1 i 
-2.31860-01 -1.8504D-0H 
-3.47340-01 -7.474B0-19i 
-6.57040-02 -8.26990-18i 
3.14350-04 +4.39590-03i 
3.14350-04 -4.39590-03i 
-5.41740-01 +0. 00000+00 i 
-5.00000-01 +0. 00000+00 i 
1.49970-03 +0.00000+00i 
-2.76000-01 +0.00000+00i 
-2.00000-01 +0.00000+00i 
-1.00000+00 +0. 00000+00 i 
-1.92000+01 +0.00000+00i 
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TABLE 2. DESIRED EIGEN VECTORS AND THEIR INTERPRETATION ”99'' STANDS FOR 
UNSPECIFIED COMPONENTS 


Measur enents 


<t 

p 

h 

h 

a 

X 

M 

/ M 

/ h 


0. 

0 . 

0 . 

0 . 

99. 

99. 

99. 

0 . 

V 

0 . 

0 . 

99. 

99. 

1. 

99. 

0 . 

99. 

Or 

0 . 

0 . 

99. 

99. 

0 . 

0 . 

0 . 

99. 

Q 

0 . 

0 . 

0 . 

1. 

99. 

0 . 

99. 

0 . 

6 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

e 

99. 

1 . 

0 . 

99. 

0 . 

0 . 

0 . 

0 . 

p 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

r 

1. 

99. 

99. 

99. 

0 . 

0 . 

0 . 

0 . 

4> 

0 . 

0 . 

i. 

99. 

99. 

0 . 

99. 

99. 

h 

0 . 

0 . 

0 . 

0 . 

99. 

1 . 

0 . 

A Thrust 

99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. - 

Actuator 

\ 

99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 

States 

99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 

/ CAS States 

99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. 


99. 

99. 

99. 

99. 

99. 

99. 

99. 

99. J 


99. 

99. 

0 . 

0 . 

99. 

99. 

1 . 

0 . /« 

99. 

99. 

99. 

99. 

0 . 

0 . 

0 . 

i. /h 
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TABLE 3. DESIRED EIGEN VALUES 


-2.0000D-01 +2-OOOOD— Oli 
-2.0000D-01 -2.0000IH)li 
-i.5000D-0i ♦i.OOOOlhOii 
-1.5000^-01 -l.OOOOD-Oli 
-i.4000IH)i +1.40000-011 
-2.50000-01 +0.00000+00i 
-1.4000H»i -1.4000D-0H 
-2.7000H)! ♦O.OOOOIHOOi 


TABLE A. OUTPUT FEEDBACK CONTROLLER 



-2.6554902 -1.7021902 -4.2989902 -1.1486902 -1.09261*02 -1.3532^*03 -1.36950*02 

1.32590-01 4.19000-41 5.W32902 2.93941*02 -0.44940*00 -4.93721*01 -7.93141*00 

-1.0067900 -5.10920*01 0 0 0 0 0 

2.10411-41 -2.31160*00 0 0 O 0 O 


-2.75031-03 


5.69630-03 

0 

0 



TABLE 5. CLOSED LOOP EIGEN VALUES 


-1.X12K02 *7.99X1-151 
-1.19651*02 *2.43381-151 
-1.00001*02 *1.11690-101 
-4.81941*41 ♦4.50581-151 
•8.80041*41 -5. 78881*41 i 
•8.80041*41 *5.78869*011 
-1.14601*41 *2.17061*011 
•1.10801*41 -2. 17061*41 i 
•4.90151*41 *1.09591-181 
-2.10241*41 *1.40319-15i 
-1.92001*41 *4.00001*40i 
-5.41871*00 *3.30731*01i 
-8.41871*40 -3.30731*011 
-5.47131*40 -4.43061*411 
-5.07131*40 *4.43061*411 


-4.13911*40 —1.40549*411 
-4.1X11*40 *1.40549*411 
-3.20281*40 *4.06319*401 
-3.20281*40 -4.06319*401 
-1.00001*00 *0.00001*40 1 
-6.94X1-01 *6.73151-011 
-1.94X1-01 -1.73151-011 
-5.00006-01 *0.00001*001 
-4.77651-01 *0.0000!*00i 
-3.40691-01 *1.15139-171 
-2.70001-01 *1.30021-171 
-2.50001-01 *1.30921-171 
-2.00451-01 -2.00379-011 
-2.00451-01 *2.00379011 
-1.5000901 -1.00001-011 
-1.5000901 41.0000MU 
-1.4000901 -1.40009-011 
-1.4000901 *1.40009011 


l desired 
' eigenvalues 

J 
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fh, M, a, n]^ 


Figure 1. Flight Test Trajectory Control System 
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MACH NO. 




Figure 3. Altitude vs. Time Response for a Ramp Mach No. Command 
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APPENDIX A- II 


AIRCRAFT - CAS LINEARIZED MODEL AND 
FTTC DESIGN USING EIGENSTRUCTURE ASSIGNMENT 
AT 1.2 MACH, 10000 ALTITUDE 


A-II-1 


1 


1 . F Matrix (31 x 31) 


LVLWHi 1 TWJ C 
7.31781' 0? -6.b€lilH0i 
-2.95580-05 -2.33420+X 
7.35430-03 -1.O45BD+02 
O.X0Q0+N O.OOOOD+X 

o.xood+x o.oxod+x 

O.XW0+X O.OXOD+X 
O.MOOD+X O.OXOD+X 
O.XXD+X C.OOXD+X 
O.XXD+X -1.29290+03 
0. WW0+X O.OXCD+X 
1.8719D-01 6.7079D+02 

O.XXD+X O.OOODIhOC 
o.ooood+x O.OXOD+X 
O.XXD+X O.OOXD+X 
-3.74370-01 -1.3416D+03 
1.87190-O2 €.?079lKDi 

1,87190+X €.70790+03 

0.00000*00 0.00000+00 
c.ooooo+oo c.ooooo+x 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0. 00000+00 
0 . 00000+00 o.ooooo+oc 
0 , 00000+00 0 . 00000+00 
0 . 00000+00 o.ooooo+oc 
0.00000+00 0.00000*00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 


COLUHK 7 THRU 12 
0 . 00000+00 0 . 00000+00 
0.00000+00 0. OOOO0+OO 

0.OO0OD+OC 0.00000+00 
0.00000+00 0.00000+00 
-9.99930-01 2.48600-02 

9.95650-01 0.00000+00 

-1.05680+00 0.00000+00 

1.19150-02 0.00000+00 

0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
1.16490+00 0.00000+00 

0 . 00000+00 0 . 00000+00 
0.00000+00 0.00000+00 
0.00060+00 0.00000+00 
-2.32980+00 0.00000+00 

M6490-O1 0. 00000+00 

1.16490+01 0.00000+00 
0 . 00000+00 0 . 00000+00 
0.00000+00 0.00000+00 
0.00000+00 C. 00000+00 
-1.05330+01 O.OOOO0+OC 
8.80000+01 0.00000+00 
0.00000+00 0.00000+00 
0 . 00000+00 0 . 00000+00 
0.00000400 0.00000+00 

0.00000+00 0.00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0.00000+00 0.00000+00 
0 . 00000+00 0 . 00000+00 
0.00000+00 0.00000+00 


-1.78931' 17 
1.00000*00 
-5.15740+00 
1.00000+00 
0.00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
O.OOOOD+OO 
3.59290+02 
-5.16670+00 
0 . 00000+00 
C. 00000+00 
-7.1B590+O2 
3.59290+01 
3.59290+03 
0 . 00000+00 
0 . 00000+00 
0.00000+00 
0.00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0.00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 


-3.71440+01 

-1.16290-06 

1.05340-03 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

O.OOOOD+OO 

1.29290+03 

O.OOOOD+OO 

2.06580- 02 
0.00000+00 
0.00000+00 
0.00000+00 
-4.13160-02 

2.06580- 03 

2.06580- 01 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
O.OOOOD+OO 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 


0. 00000+ 00 
0 . 00000*00 
0 . 00000+00 
O.OOOOD+OO 
-5.07930-01 
-1.13270+02 
3.00490+01 
0.00000+00 
0.00X0+00 
O.OOOCD+OC 
-1.32530+02 
0.00X0+00 
0.00X0+00 
0.00000+00 
2.65050+02 
-1.32530+01 
-1.32530+03 
O.OOOCD+OC 
0.00000+00 
0.00000+00 
6.25M0+O1 
0.00000+00 
0.00000+00 
O.OOOOD+OO 
O.OOOOD+OO 

0.00000+ oc 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 


O.OOOOD+OO 
0.00000+00 
0.00000+00 
O.OOOOD+OO 
1.19140-02 
-5. 06530+00 
-6.0755H3 
1 . 00000+00 
0. 00000+ OC 
O.OOOOD+OO 
-5 92640+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOCD+OC 
1. 18530+01 
-5.9264D-01 
-5.92640+01 
O.OOOOD+OO 
-3.56130+02 
0.00000+00 
2.66510+00 
O.OOOOD+OO 
0. OOOOD+OO 
O.OOOOD+OO 
O.XXD+X 
O.OOOOD+X 
O.OOXD+OO 
O.XX0+OO 
0.00X0+00 
O.OOOCD+X 
O.OOXD+OC 


1.54400-03 

0.0X00+00 

-5.03520-08 

O.OXO0+X 

0.0X00+00 

O.OXO0+X 

0.0X00+00 

0.0X00+00 

O.OXO0+X 

0.00X0+00 

-1.03200-06 

O.OXO0+X 

O.OXO0+X 

O.OXO0+X 

2.06400-06 

-1.03200-07 

-1.03200-05 

O.OXO0+X 

O.OOOOD+OO 

O.OOOOD+X 

O.OOOOD+X 

O.OOOOD+X 

O.OXO0+X 

O.OOOOD+X 

O.OXO0+X 

O.OXO0+X 

O.OOOOD+X 

O.OXO0+X 

O.OOXD+X 

C.OXO0+X 

O.OOOOD+X 


3.79820-01 
-3.48860-06 
2.91380-06 
O.OXO0+X 
O.OXOD+X 
O.OOXD+X 
O.OXO0+X 
O.OOXD+X 
O.OOOOD+X 
-2.00X0-01 
-5.36670-05 
O.OOXD+X 
O.OXO0+X 
O.OXOD+X 
1.07330-04 
-5. 36670-06 
-5.36670-04 
O.OXO0+X 
O.OOXD+X 
O.OOXD+X 
O.OOXD+X 
O.OOXD+X 
O.OXOD+X 
O.OOXD+X 
O.OXOD+X 
O.OOXD+X 
O.OOXD+X 
O.OXOD+X 
O.OOXD+X 
O.OXOD+X 
O.OOXD+X 


O.OOXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 

O.OOXD+X 

O.XXD+X 

-5.XXD+01 

O.OOXD+X 

1.50XD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 

O.OOXD+X 

O.XXD+X 

O.XXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 


O.XXD+X 

O.XX0+X 

O.XXD+X 

O.XXD+X 

O.XXD+X 

O.XXD+X 

O.XXD+X 

O.XXD+X 

O.XXD+X 

O.XXD+X 

3.00XD+01 

-3.33330-01 

O.XXD+X 

O.XXD+X 

-6.XXD+01 

3.00XD+X 

3.XXD+02 

O.OOXD+X 

O.XXD+X 

O.OOX0+X 

O.OOXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.OOXD+X 

O.XXD+X 

O.OOXD+X 
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F Matrix (Cont'd) 


COt URK; 13 T*D 18 
O.OOW5+W 0.00005+00 

0.00005+00 o.oooor»w 

0.00005+00 0.00005*00 

0.00005*00 0.00005*00 

0.00001*00 0. 00000*00 
O.W0O5*W O.OW05+W 
0.00000*00 O.OWDIhW 
0.W005+W O.OOOONOC 
0.00000*00 0.00005*00 

0.00 OOP* 00 0.00000+00 

0.WW5+W O.OW05+W 
0.00000*00 O.OW05+OO 
0.00005*00 0.00001*00 

1,00000*02 -1.00000*0? 
e.oooop*oo o.ooooihoc 
O.OOOOP+OO 0.00005*00 
O.W0O5+W 0.00000*00 
0.00000*00 3. 49045-01 

O.OOOOIHOC 0.00000+00 

0.00000*00 o.ooooc^oo 

O.OOOOP+OO C.OOOOP+OO 
0.00005+W C. 00000+00 
0.00000*00 O.OW05+W 
0. 00000+00 0.00005+00 

0.00000+00 O.OOOOIHOC 
C.W005+W 0.00005*00 

0.00000*00 O.OOOOIHOC 
0.00005+W O.OOOOIHOC 
0.00005*00 0.00005*00 

0.00005*00 0.00005+00 

0.00005*00 0. 00005+00 


COLUHKS 19 1WIU ?4 
0.00005*00 0.00005+00 

0.00005+M 0.00005+00 

0.00005+00 0.00005+00 

0.00005+00 0.00005+00 

0.00005+00 -2.088 Sr-02 
0.00005+00 5.04585+01 

0.00005+00 3.73455+00 

0.00005*00 0.00005*00 

0.00005*00 O.OOOOIHOC 
0.00005+00 0.00005*00 

0.00005*00 6. 60565+01 

0.00005*00 0.00005+00 

0.00005+00 0.00005+00 

0.00005*00 0.00005*00 

0.00005+00 -1.32115+02 
0.00005+00 €.€0555+00 

0.00005+00 €.€0565+02 

0.00005+00 0.00005+00 

•1.00005+02 O.OOOOIHOC 
€.96085-01 -2.00005+01 
0.00005+00 -€.€3995+00 
0.0000 5+ OC 0.00005+00 

0.00005*00 0.00005*00 

0.00005*00 0.00005*00 

0.00005*00 0.00005+00 

0.00005+ OC 0.00005+00 

0.00005+00 0.00005+00 

0.00005+00 0.00005+00 

C. 00005+00 0. 00005+00 

0.00005+00 0.00005+00 

0.00005+00 0.00005+00 


0.00005+00 
0.00005+00 
0.00005*00 
0.00005+00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
1.00005*01 
0.00005*00 
0.00005*00 
0.00005*00 
-3.00005*01 
1.00005*00 
1.00005*02 
0.00005+00 
0. 00005+00 
0. 00005+00 
O.OOOOD+OO 
0.00005+00 
0.00005*00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005*00 

-€.€0005+01 

0.00005+00 

0.00005*00 

0.00005*00 

3.43805*03 

0.00005*00 

0.00005*00 

0.00005*00 

0.00005*00 

0.00005*00 

0.00005*00 


0.00005*00 

0.00005+00 

0.00005*00 

0.00005+00 

0.00005+00 

0.00005*00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

1.00005+01 

0.00005*00 

0.00005+00 

0.00005+00 

O.OOOOD+OO 

0.00005+00 

1.00005+02 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005*00 

0.00005*00 


0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005+0: 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005*00 
0.00005*00 
0.00005*00 
-8.80005*01 
-5. 00005-01 
-5.00005-01 
3.43805+03 
0.00005+00 
0.00005*00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


0.00005*00 
0.00005+00 
0. 00005+00 
0.00005+00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.0CW5+W 
0.00005+00 
0.00005*00 
0.00005*00 
0.00005+00 
0.00005*00 
-1.00005+02 
3.49045-01 
0.00005+00 
0.00005*00 
0.00005+00 
0. WOO 5* 00 
0.00005*00 
0.00005*00 
0.00005+00 
C. WOO 5* 00 
0.00005*00 
0. WOO 5+ 00 
O.OOW5+W 
0.00005+00 
0.00005*00 


0.00005*00 
0.00W5+W 
0.00005+W 
0. 00005* W 
0.00W5+W 
0.00005+W 
O.WD05+W 
0.W005+W 
0.WW5+W 
O.OOW5+W 
0.WW5+W 
O.OOW5+W 
0.0CW5+W 
0.00005+W 
O.OOW5+W 
0.00005+W 
0.00005+W 
0.WW5+W 
0.00005+W 
0.00W5+W 
O.OOW5+W 
O.OOW5+W 
-5.00005-01 
-5.00005-01 
3.43805*03 
0.00W5+W 
O.OOW5+W 
0.00W5+W 
O.OOW5+W 
0.00005+W 
0.00005+W 


-9.01B45+W 
-2.38305-01 
-4.51535*01 
0.00005+W 
0.WW5+W 
0.W005+W 
O.WD05+W 
O.W0O5+W 
0.00005+W 
0.WDO5+W 
-6.38105+02 
0. W0C5+W 
0.W005+W 
0.WW5+W 
1.27625+03 
-6.38105+01 
-6.38105+03 
-2. WOO 5+01 
0.WD05+W 
0.00005+W 
0.W0O5+W 
O.W005+W 
0.W005+W 
0.W005+W 
0.00005+W 
O.W005+W 
C.W005+W 
0.00005+W 
0.00005+W 
0.W005+W 
O.W0O5+W 


O.OOW5+W 

0.WW5+W 

C.W0C5+W 

0.WW5+W 

O.W0O5+OO 

0.WW5+W 

0.W005+W 

0.00005+W 

0.00W5+W 

0.W0O5+W 

0.00005+W 

0.00W5+W 

0.00005+00 

0.00005+00 

0.00005*00 

0.00005+00 

0.00005+00 

0.00005*00 

0.00005+00 

0.00X5+00 

0.00005+00 

0.00005*00 

0.00005*00 

-5.00005-01 

3.43805*03 

O.WW5+OC 

0.00005*00 

0.00005*00 

0.00005*00 

0.00005*00 

0.00005*00 


A-II-3 


F Matri x (Cont 1 d) 


CDIUKKS » WU 30 

o.ooood+oc c.ooood+oc 
0.00000+00 0.00000+ 00 
0.00000+00 O.OOOOD+OO 
0.00000+00 O.OOOOD+OO 
0.00000*00* O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO C.OOOOD+OO 
O.OOOCD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OC 
O.OOOOD+OC O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OC 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO C.OOOOD+OO 
O.OOOOD+OC O.OOOOD+OC 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
"0.00000+01 O.OOOOD+OO 
5.0000D-01 O.OOOOD+OO 
5.5000D+01 5.S000D+01 
O.OOOOD+OC O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OC O.OOOOD+OC 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

C.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

C.OOOOD+OO 

C.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

C.OOOOD+OO 

C.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

C.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-5.5000P+C1 

l.OOOOD+O? 

O.OOOOD+OO 

C.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-l.OOOOD+O? 

4.B866D-01 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-3.8962D-02 

-3.9064D+0C 

-U634D+01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-4.5706D+00 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

9.1410D+00 

-4.57051-01 

-4.5705D+01 

O.OOOOD+OO 

O.OOOOD+OC 

O.OOOOD+OC 

-1.2K3M? 

O.OOOOD+OC 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OC 

O.OOOOD+OC 

O.OOOOD+OO 

O.OOOOD+OC 

-2.B000D+01 

O.OOOOD+OC 

O.OOOOD+OO 


O.OOOOD+OC 

O.OOOCD+OO 

C.OCOCD+OO 

O.OCOCD+OC 

O.OOOOD+OC 

O.OOOOD+OO 

O.OOOOD+OC 

O.OOOOD+OC 

O.OOOOD+OC 

O.OOOOD+OC 

O.OOOOD+OC 

O.OOOOD+OC 

O.OOOOD+OC' 

C.OOOCD+OC 

O.OOOOD+OC 

C.OOOCD+OC' 

O.OOOOD+OC 

O.OOOOD+OC 

C.OOOCD+OC 

O.OOOOD+OO 

O.OOOOD+OC 

C.OOOCD+OC 

O.OOOCD+OO 

O.OOOOD+OC 

C.OOOOD+OC 

C.OOOOD+OO 

C.OOOCD+OC 

O.OOOOD+OC 

O.OOOOD+OC 

-l.OOOOD+OC 

I.S200D+01 


COLUHH 31 
O.OOOOD+OC 
O.OOOOD+OO 
0. 00000+00 
O.OOOOD+OC 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OC 
O.OOOOD+OC 
O.OOOOD+OC' 
O.OOOOD+OO 
O.OOOCD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OC 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
-5.9000D+01 
O.OOOOD+OO 
O.OOOCD+OO 
O.OOOOD+OC 
-1. 52000+01 


A-II-4 


2 


G-Matrix (31 x 4). 


O.OOOOD+OO 

O.OODOIHOO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOCD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

2.0000D-01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

l.OOOOD+Ol 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

i.OOOOD+OO 

1.0000D+02 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

1.0000D+02 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

3.4380D+03 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


A-II-5 


3. H-Matrix (1? x 3D 


COLUMNS 
0. OOOOD+OO 
1.2151D-03 
0. OOOOD+OO 
7. 3543D-03 
0. OOOOD+OO 
O.OOOOD+OO 
0. OOOOD+OO 
0.00000+00 
0.00000+00 
9.28150-04 
0.00000+00 
1.18920-03 
0.00000+00 
0.00000+00 
0.00000+00 
9.99540-01 
-2.26010-03 


COLUMNS 

9.95650-01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-1.05680+00 

1.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


1 THRU 6 
0. OOOOD+OO 
9.3810D+01 
0. OOOOD+OO 
-1.0458D+02 
0. OOOOD+OO 
0. OOOOD+OO 
0. 00005+00 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
1. OOOOD+OO 
9.2813D+01 
0. OOOOD+OO 
0. OOOOD+OO 
-1.2929D+03 
-1.5381D+01 
-6.1083D-01 


7 THRU 12 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 
O.OOOOD+OO 
0. OOOOD+OO 
1. OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


O.OOOOD+OO 
O.OOOOD+OO 
1. OOOOD+OO 
-5.1574D+O0 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-5.0353D-08 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

9.7656D-01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

4.7963D-06 


O.OOOOD+OO 
-3.1062D-05 
O.OOOOD+OO 
1.0534D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
-1.2611D-06 
O.OOOOD+OO 
1. OOOOD+OO 
1.2929D+03 
O.OOOOD+OO 
8.8325D-06 


O.OOOOD+OO 

-3.7795D-06 

O.OOOOD+OO 

2.9137D-06 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-3.7834D-06 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

1.1805D-02 


-1.13270+02 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

-2.0410D+01 

3.0049D+01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


-5.0653D+00 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
1. OOOOD+OO 
O.OOOOD+OO 
-6.0755D-03 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 
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H-Matrix (Cont'd) 


COLUMNS 13 THRU 18 
O.OOOOD+OO 0.0000D+00 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 


CXUHNS 19 THRU 24 
O.OOOOD+OO 5.6458D+01 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO -8.3936D-01 
O.OOOOD+OO 3.7345D+00 
O.OOO.D+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 
9.5783D+00 
O.OOOOD+OO 
-4.5153D+01 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
9.5847D+00 
O.OOOOD+OO 
OjOOOOD+OO 
O.OOOOD+OO 
O.OOOOD+OO 
-1. €€160-01 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 
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H-Matrix (Cont'd) 


COLUHKS 25 TWU 30 
0.00000+00 0.00000+00 
0.00000+00 0.0000D+00 

0. 00000+00 0.0000D+00 

0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 


0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


0 . 00000+00 - 3 . 00640+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0 . 00000+00 - 1 . 56560+00 
0 . 00000+00 - i . 16340+01 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 
0 . 00000+00 0 . 00000+00 


0.00000+00 

0.00000+00 


0.00000+00 

0.00000+00 


0.00000+00 

0.00000+00 


0.00000+00 

0.00000+00 


0.00000+00 

0.00000+00 


COL WIN 31 

0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 


0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000+00 

0.00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 * 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 
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Controller Design 


U = GAIN * z 


where 


z are perturbed values of 
[<f> p h h ax M /M /h]^ 


and 


u ■= f 6T Se 6a 6Y ] T 
ap ap ap 


CAIN 


COLUMNS 
-2.65549-02 
1.32599-01 
-6. 00679+00 
2.60419-01 


l mu 6 
-6.70219-02 
4.19069-01 
-5.10929+01 
-2.31169+00 


-4.29899-02 

5.96329-02 

3.71649-12 

4.06779-13 


-1.14889+02 

2.93949+02 

1.76329-08 

1.93219-09 


-1.09269+02 

-6.44949+00 

-7.08129-12 

-4.52929-13 


COLUMNS 7 Wii 8 
-1.38959+02 -2.75039-03 
-7.93149+00 5.89639-03 

-1.83909-11 3.67529-13 

-1.01629-12 4.13119-14 


-1.35329+03 

-4.93729+01 

€.16619-10 

7.29119-11 
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Open-loop ei g env al ues 


-5.94560‘01 +2.02930+011 
-5.94560+01 -2.02S30+OU 
-1.36740+0? -3.2679D-15) 
-4.86730+01 *3.O4340-17i 
-6.77^8IKC1 +5. 82180+01 I 
-8.77580+01 -5.82180+011 
-J.C1M0+02 +5.17230-151 
-1 .OOOOIhO? +4.21440-161 
-8.819SW1 +3.54E9H6i 
-2.32P60+01 +8.24J60-17i 
-7.41170+00 +3.35740+OH 
-7.411717*00 -3.35740+Oli 
-7.97251*00 +8.2€550+00i 
-7.97250+00 -6. 20550+00 i 
-1.21740+01 +1.01920+01i 
-1.21760+01 -1.01920+011 
-2.94380+00 +5.31130+001 
-2.94380+00 -5.31130+001 
-2.31960-01 +1.85040-011 
-2. 31961-01 -1.86040-011 
-3.47340-01 -7.47480-191 
-6.57040-02 -8.26990-181 
3.14350-04 +4.39590-031 
3.14350-04 -4.3!590-03i 
-6.41740-01 +0.00000+001 
-5.00000-01 +0.00000+001 
1.49970-03 ♦0.00000*001 
-2.76000-01 +0.00000+001 
-2.00000-01 +0.00000*001 
•1.00000+00 +0.00000*001 
-1.92000+01 +0.00000+001 


Closed-loop eigenvalues: 


-1.36120+02 +7.99360-15i 
-1.19550+02 +2.43380-151 
-1.00000+02 +6.11890-161 
-8.81940+01 +4.50580-151 
-8.80040+01 -5.78860+011 
-8.80040+01 +5.78860+flli 
-6.10800+01 +2.17060+0H 
-6.10600+01 -2.17060+011 
-4.90650+01 +1.09590-16i 
-2.60240+01 +6.4C3tO-15i 
-1.92000+01 +0.00000+00i 
-5.41870+00 +3. 30730+01 i 
-6.41870+00 -3. 30730+01 i 
-5.07130+00 -4.43060+01i 
-5.07130+00 +4. 43060+011 
-4.13910+00 -1.40540+011 
-4.13910+00 +1.40540+011 
-3.20280+00 +4.06310+00i 
-3.20280+00 -4.06310+00i 
•1.00000+00 +0.00000+001 
-6.94390-01 +6.73150-011 
-6.94390-01 -6.73150-Cli 
-5.00000-01 +0.00000+001 
-4.77650-01 +0.00000+001 
-3.40690-01 +1.15130-171 
-2.70000-01 +1.30020-17i 
-2.50000-01 +1.30920-171 
-2.00450-01 -2.00370-011 
-2.00450-01 +2.00370-011 
-1.50000-01 -1. 00000-01 i 
-1.50000-01 +1.00000-011 
-1.40000-01 -1. 40000-01 i 
-t. 40000-01 +1.40000-011 
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CLOSED-LOOP TIME RESPONSE 
MACH NO VS TIME FOR A RAMP MACH NO COMMAND 
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ALTITUDE 



TIME 


CLOSED-LOOP 
ALTITUDE VS TIME FOR 


TIME RESPONSE 
A RAMP MACH NO COMMAND 
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APPENDIX A-III 


AIRCRAFT - CAS LINEARIZED AND 
FTTC DESIGN USING EIGENSTRUCTURE ASSIGNMENT 
AT 0.7 MACH, 10000 ALTITUDE 


1 . F MATRIX (31 x 31) 


cauws 
1.177BM? 
-1.24925 04 
7.6?7ll*04 
0.00005+00 
5.34815-13 
-3.47925-09 
-1.3023M0 
0.00001*00 
0.00005*00 
0.00001*00 
1.0380001 
0.00001*00 
0.00000+00 
0.000004 00 
-2.C759M1 
1.03800-0? 
1.03801*00 
0.00001*00 
0.00001*00 
0.00001*00 
1.1261D-09 
0.00001*00 
0.00000400 
0.00005+00 
0.00000+00 
C.00005+00 
0.00005+00 
0.00005+00 
O.OOttO+OO 
0.00005+00 
0.00005+00 


1 THRU 0 
2.21255+01 
-1. 32925* 00 
-6. €2345*00 
0.00005+00 
71C7D 10 

B. 9156MQ 
-4.52085-08 

0.00005+00 

-7.54185*02 

0.00005+00 

7.97755+02 

0.00005+00 

0.00005+00 

C. 00005+00 
-1.55555+03 

7.97755+01 
7.97755+03 
0. 00005+00 
0.00005+00 
0.00005+00 
-5.54165-07 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0. 00005+00 
0.00005+00 


O.OOOCD+OO 
1.00005+00 
-3.12485+00 
1.00005+00 
O.OOOOD+OO 
0.00005+00 
0.00005*00 
0.00005+00 
0.00001*00 
0.00005+00 
4.00905+0? 
-5.16675*00 
0.00005+00 
0.00005+00 
-6.01915+02 
4.00965+01 
4.00965+03 
0.00005+00 
0.00005+00 
C. 00005+00 
0.00005+00 
0.00005+00 
0. 00035+00 
0.00005+00 
0.00005+00 
C. 00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


-3. 21445+01 
0.00005+00 
1.09765-03 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
C. 00005+00 
7. 54 181*02 
0.00005+00 

2.24955- 0? 
C. 00005+00 
0.00005+00 
0.00005+00 

-4.49915-02 

2.24951*03 

2.24955- 01 
0.00005+00 
0.00005+00 
0.00005+00 
C. 00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
C. 00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
-2.88985-01 
-4.3656 5+01 
1.05735+01 
0.00005+00 
0. 00005+00 
0.00005+00 
-5.10785+01 
0.00005+00 
C. 00005+00 
0.00005+00 
1.02165+02 
-5.10785+00 
-5. 10785+02 
0.00005+00 
C. 0000 5+00 
0.00005+00 
2.97005+01 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0. 00005+00 
0.00001*00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


cauws 

0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
-9.99445-01 
1. 54475+00 
-7.50535-01 
3.34625-02 
0.00005+00 
0.00005+00 
1.80735+ 00 
0.00005+00 
0.00005+00 
0.00005+00 
-3.61465+00 
1.80735-01 
1.80735+01 
0.00005+00 
0.00005+00 
C. 00005+00 
-6.00595+00 
6.80005+01 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


7 TWU 12 
0.00005+00 
O.OOOOD+OO 
0.00005+00 
0.00005+00 
4.25975-02 
-2.56495-23 
3.23125-24 
0.00005+00 
0.00005+00 
0.00005+00 
2.56491*23 
0.00005+00 
0.00005+00 
0.00005+00 
-5,16995-23 
0.00005+00 
-2.06805-2 2 
0.00005+00 
0.00005+00 
0.00005+00 
-6.05945-07 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


0.00005+00 
O.OO0O5+0O 
C. 00«5*00 
0.00005+00 
0.00005+00 
0.00005+00 
O.OOOOD+OO 
0.00005+00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00001*00 
0.00005*00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005+00 
0.00005+00 
0.00005+00 


3.79245-01 

-1.61575-05 

6.09625-06 

0.00005+00 

0.00005+00 

0.00005*00 

0.00005+00 

0.00005+00 

0.00005*00 

-2.00005-01 

-3.04265-05 

0.00005+00 

O.OOOOD+OO 

0.00005+00 

6.00535-05 

-3.04265-06 

-3.04261*04 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

6.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 


0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005+00 

0.00005*00 

0.00005+00 

-5.00005+01 

0. 00005+00 

1. bOOOD+OC 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 


0.00005+00 
0.00005+00 
0.00005*00 
0.00005*00 
3. 34435-02 
-2.62555*00 
-4 . 85921-02 
1.00005+00 
0.00005*00 
0.00005*00 
-3.30585*00 
0.00005*00 
0.00005*00 
0.00005+00 
6.61165*00 
-3.30585-01 
-3.30505*01 
0.00005*00 
-3.58135*02 
0.00005*00 
1.05905*00 
0.00005*00 
0.00005+00 
0.00005*00 
0,00005+00 
0.00005*00 
0.00005*00 
0.00005*00 
0.00005*00 
C. 00005*00 
0.00005*00 


0.00005+00 
0.00005+00 
0.00005*00 
0.00005+00 
0.00005+00 
0.00005*00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
3.00005+01 
-3.33335^01 
0.00005+00 
0.00005+00 
-6.00005+01 
3.00005+00 
3.00005+02 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
C. 00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
0.00005+00 
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F Matrix (Cont'd) 


COtUWS 13 THRU 18 
O . OOOOO+OC 0.00000+00 
0.00000*00 0.00000*00 
C . OOOOO+OO o . ooooo+oe 
O . OOC * OIH ( K > 0.00001*00 

0.00000*00 o . ooooo+oo 
0. 00000+00 0.00000+00 
0,00000*00 0.00000+00 
0.00000+00 C . 00000+00 

o.ooooiwoo o . ooooo*oo 
0 . 00001*00 o . oooot+oo 
0.00000+00 0.00001*00 
0.00000+00 0.00001*00 
0.00001*00 0.00000+00 
1.00000+0? -1.00001*02 
0.00001*00 0.00001*00 
0.00000+00 0.00000+00 
0.00001*00 0.00001*00 
0.00000+00 3.4*041*01 

o . oowo*oo 0 . 00000+00 
0.00001*00 0.00001*00 
0.00000+00 0.00000+00 
0.00000+00 C . 00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+ 00 0.00000+00 
0. 00000+00 C. 00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 


C 0 LUW 5 19 THRU ?4 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 -2.36270-02 
0.00000+00 1.90290+01 

0.00000+00 2.11050+00 

0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 c . ooooo+oo 

0.00000+00 2.22630+01 

0.00000+00 0. 00000+00 
0.00000+00 0.00001*00 
0.00000+00 0.00000+00 
0.00000+00 -4.45270+01 
0.00000+00 2.22630+00 

0.00000+00 2.22630+02 

0.00000+00 0.00000+00 
-1.00000+02 0.00000+00 
6.9B06D-01 -2.00000+01 
0.00000+00 2.06901*00 

0.00001*00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00001*00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000*00 
0.00000+00 0.00000+00 


O.OOOOPtOC 
0 . 00000*00 
0.00001*00 
0 . 00000*00 
0 . 00000*00 
0 . 00000*00 
0 . 00000*00 
0 . 00000*00 
0 . 00000*00 
0.00001*00 
1 . 00000*01 
0.00000+00 
0 . 00000*00 
0 . 00000*00 
-3.00000+01 
1.00000+00 
1.00000+02 
0.00000+00 
0.00000+00 
c . 00000+00 
0.00000+00 
c . ooooo+oo 
c . ooooo+oo 
c . ooooo+oo 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0 . 00000*00 
0.00000+00 
0.00000*00 


0.00000+00 
0.00000+00 
0.00000+00 
0.00001*00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00001*00 
0.00000+00 
0.00000+00 
0.00001*00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
-€.60000+01 
0.00001*00 
0.00000+00 
0.00000+00 
3.43800+03 
0.00000+00 
0.00001*00 
0,00001*00 
0 00000+00 
0.00000+00 
0.00000+00 


0 . 00000*00 
0.00000*00 
c . ooooo+oo 
0.00000*00 
0. 00000+00 
0.00000+00 
0 . 00000*00 
0. 00000+00 
c . ooooo+oo 
c . ooooo+oo 
1.00000*01 
0.00001*00 
0.00000+00 
0.00001*00 
0 . 00000*00 
0.00000+00 
1.00000+02 
C. 00000+00 
C. 00000+00 

c . ooooo+oo 

c . ooooo+oo 

0.00000+00 

0.00000*00 

0.00000*00 

0.00001*00 

c . ooooo+oo 

0 . 00000*00 

c . ooooo+oo 

0 . 00000*00 

0.00000*00 

0.00000*00 


c . ooooo+oo 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000*00 

0.00000*00 

0.00000+00 

0.00000+00 

c . ooooo+oo 

0.00000+00 

O . OOODO+OO 

c . ooooo+oo 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-6.00000+01 

-6.00000-01 

-5.00000-01 

3.43800+03 

0.00000+00 

0 . 00000*00 

0 . 00000*00 

o . ooooo+oo 

0 . 00000*00 

0.00000+00 


0. 00000+00 
0.00000*00 
0.00000+00 
0.00000+00 
0.0000 0+00 
0.00000*00 
0.00000+00 
0. 00000+00 
0.00000*00 
o.oooco+oc 
o.ooooo+oc 
0.00000+00 
0.00000+00 
0.00000+00 
0.00001*00 
0.00000+00 
-1.00000+02 
3.49040-01 

o . ooooo*oc 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 

0.00000*00 


0.00000*00 

o . ooooo+oc 

0.00000+00 

c . ooooo+oo 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000*00 

o . ooooo+oc 

c . ooooo+oo 

0 . 00001*00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-5.00000-01 

-5.00000-01 

3.43800+03 

0.00000+00 

c . ooooo+oo 

c . ooooo+oo 

o . ooooo+oc 

o . ooooo+oc 

0.00000+00 


-1.40520+00 

-1. B 8490-01 

-1.94730*01 

0.00000*00 

0.00000*00 

0 . 00000*00 

0 . 00000*00 

0.00000+00 

0.00000*00 

0 . 00000*00 

-2.66800*02 

0 . 00000*00 

0 . 00000*00 

0.00000+00 

5.33190*02 

-2.66600+01 

-2.66600+03 

-2.00000*01 

0.00000*00 

0 . 00000*00 

0 . 00000*00 

0 . 00000*00 

0.00000+00 

o . ooooo+oc 

0.00000*00 

0.00000+00 

0.00001*00 

0.00000*00 

c.ooooo+oo 

0.00000+00 

6.00000*00 


6.00000+00 

o . ooooo+oc 

6.00000+00 

c . ooooo+oo 

6.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

o . ooooo+oc 

0.00000+00 

0.00000+00 

0.00001*00 

o . ooooo+oc 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

c . ooooo+oo 

0.00000+00 

0.00000+00 

-5.00000-01 

3.43800+03 

0.00000+00 

o . ooooo+oc 

0.00000+00 

0.00000+00 

o . ooooo+oc 

0.00000+00 
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F Matrix (Cont'd) 


COLWWS 25 THRU 30 
O . OOOOD+OQ 0. 00000*00 

0.00000+00 

0. 00000+00 

0.00000+00 

0.00000+00 

COLUMN 31 
0.00000+00 

0.00000+00 

0.00000+00 

0. 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0. 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0. 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

C . 00000+00 

0.00000+00 

0.00000+00 

0. 00000+00 

-5.10770-02 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-2.13100+00 

o . ooooo+oo 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-5.14510+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

O . OOOOD+OO 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-2.49330+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

4.98660+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-2.49330-01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-2.49330+01 

0.00000+00 

o .- ooooo+oo 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0. 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-6.41050+01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-€.00000+01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

5.00000-01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

6.90000+01 

5.90000+01 

-5.90000+01 

0.00000+00 

0.00000+00 

0.00000+00 

-5.90000+01 

0.00000+00 

0.00000+00 

1.00000+02 

-1.00000+02 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

4 . 88€€ D— Oi 

-2.80000+01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-1.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

1.92000+01 

-1.92000+01 
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2 


G-Matrix (31 x *0 


O.OOOOIHOO 

0.00000+00 

0.0000D+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

O.OOOOD+OO 

O.OOOOIHOO 

2.00000-01 

0.00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

1.00000+01 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

1 . 00000+00 

1 . 00000+02 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 


0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

1.00000+02 

0.00000+00 

0.00000+00 

O.OOOOIHOO 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


0.00000+00 

O.OOOOIHOO 

0.00000+00 

0.00000+00 

0.00000+00 

O.OOOOIHOO 

0.00000+00 

O.OOOOIHOO 

O.OOOOIHOO 

O.OOOOIHOO 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

3.43800+03 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 
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I 


3. H-Matrix (17 x 31) 


COLUMNS 1 THRU C 
-3.4792 D -09 9.09490-10 

2.93880-03 3.11170+01 

O . OOOOD+OO 0.00000+00 
7.62710-04 -6.62340+00 
0.00000+00 0.00000+00 
3.52460-11 -6.36820-09 
-1.30230-10 -4.92070-08 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
9.28150-04 0.00000+00 

0.00000+00 1.00000+00 
2.93090-03 3.10420+01 

0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 -7.54180+02 
9.99430-01 -2.52080+01 
-2.67930-04 1.72930+00 


COLUMNS 7 THRU 12 
1.54470+00 0.00000+00 

0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 -5.82080-08 
-7.58530-01 0.00000+00 

1.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 1.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 
0.00000+00 0.00000+00 


0.00000+00 

5.89490-05 

1.00000+00 

-3.12480+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

2.91470-05 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-2.14200-08 


0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000+00 

0.00000+00 

0 . 00000+00 

9.76560-01 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 


0.00000+00 

0.00000+00 

0.00000+00 

1.09760-03 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

1.00000+00 

7.54180+02 

0.00000+00 

5.58790-06 


0.00000+00 

-5.17910-06 

0.00000+00 

6.09620-06 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-5.18630-06 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

1.17920-02 


-4.36560+01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

-6.77380+00 

1.05730+01 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0. 00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 


-2. 82550+00 
0.00000+00 
0.00000+00 
0.00000+00 
1.00000+00 
0.00000+00 
-4.85920-02 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 


0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 
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H-Matrix (Cor'ii'd) 


OIUWS 13 THRU 18 

o.ooood+oo o.oooor^oo 

O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
0.0000 IK 00 O.OOOOD+OO 

O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 


CQLintS 19 THRU 24 
O.OOOOD+OO 1.9029D+01 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO -5.5384D-01 
O.OOOOD+OO 2.1105D+00 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 
O.OOOOD+OO O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 


O.OOOOD+OO 

4.4173D+00 

O.OOOOD+OO 

-1.9473D+01 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

4.4225D+00 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

1.0411D-01 


O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 

O.OOOOD+OO 
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H-Matrix (Cont'd) 


caums 25 THRU 30 
0.00000+00 O.OOOOD+OO 

o.ooooo+oo o.ooooo+oo 
0 . 00000+00 


0 . 00000+00 

O.OOOOIKOO 

0 . 00000+00 

Q. 00000+00 

o.ooooo+oo 
O.OOOOIKOO 
0. OOOOO+OQ 
0 . 00009+00 
0.00009+00 
0.00009+00 
0.00009+00 
t.MOOMO 
0.00009+00 

o.ooooo+oo 

0.00009+00 


0 . 00000+00 

0 . 00000+00 

o.ooooo+oo 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 


0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


-2.13100+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

-1.19730+00 

-5.14510+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 


can* 3 i 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0 . 00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 
0.00000+00 


0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0 . 00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 

0.00000+00 
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Con tro ller Des i gn 


U - Gai n*z 


where 


z are perturbed values of 
[i p h h H /H / h] T 

and 


[ 6T 6e 6a 6Y ] 
ap ap ap 


C0LUHK5 1 THRU 6 
-1.47779-02 -7.30969-02 
1.19249-02 3.29349-02 

-2.54379+00 -1.29559+01 
1.41039-01 -5.31879-02 


CAIK 


COLUKS 

-1.39169+02 

-2.47919+00 

-1.37029-06 

-2.37389-06 


7 THRU 8 
-7.26979-03 
1.07329-03 
3.71979-12 
-1.14939-13 


-9.28509-02 

1.01559-02 

2.92559-11 

-5.86839-12 


-1.71419-01 

3.91519-02 

4.36199-11 

2.93279-11 


-1.35599+02 

-2.04429+00 

7.25459-08 

-5.98109-08 


-1.53529+03 

-1.96199+01 

-5.46979-06 

-7.43389-07 
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Open- l oop ei genval ues 


-1.28843+0? -8.88i8M6i 
-8.54286+01 -4.33bl6+0h 
-8.54206+01 +4. 93516+01) 
-8.80396+01 +5.38190-l?i 
-1.00600+0? + 1.1917M8) 
-5.02626+01 -1.11456-i5i 
-0.81750+01 +1.88930-101 
-1.00006+0? +?.B7306-15i 
-2.88876+01 -3.60876+00) 
-2.88S76+01 +3. 60970+001 
-1.46890+01 -2.78070+01) 
-1.46890+01 +2. 78070+01 i 
-1.53710+01 -4. 01640-16i 
-3.83896+00 -2.47B10+00) 
-3.83896+00 +2.47816+00i 
-7.83276+00 +1.05590-16i 
-1.77106+0© +3.14536+00) 
-1.77100+00 -3.14530+0©i 
-4.88230-01 +7.39970-18) 
-8.12466-01 +1.71020-17) 
-2.74380-01 47.87780-021 
-2,74380-01 -7.87780-021 
-6.00000-01 -2.3109M7) 
-1.21340-02 +1.447B0-16i 
-5.22000-11 45.80230-151 
1.82080-03 -5.8158M51 
7.85730-03 +1.40940-17i 
-1.71390-01 -5.18270-18) 
-2.00000-01 40.00000+001 
-1.00000+00 40.00000+00) 
-1.32000+01 40.00000+00) 


Closed-loop eigenvalues 


-1.28830+02 +0.00000+001 
-1.02840+02 -1.840l*hl5i 
-1.00003+02 -4.2141M0i 
-8.81740+01 -2.8I84M51 
-8.78710+01 -2.B9C26-151 
-8.54410+01 44.33370+Oli 
-8.54410+01 -4. 33370+01 i 
-5.02700+01 48.03933-18) 
-2.83840+01 -3. 84800+00 i 
- 2 . 83843+01 + 3 . 64600 + 00 ) 
-1.32000+01 +0.00000+00i 
-1.47530+01 -2. 78530+01 i 
-1.47590+01 +2.78530+01i 
-1.00670+01 +1.27500+01) 
-1.00670+01 -1.27500+Oli 
-3.87210+00 +1.87720+00) 
-3.87210+00 -1.87720+OOi 
-2.10110+00 +3.18380+001 
-2.10110+00 -3.16380+OOi 
- 1 . 00000+00 + 0 . 00000 + 00 ) 
-5.54010-01 +0.00000+00) 
-5.00000-01 +0.00000+00i 
-3.90830-01 -1.3330H7i 
-2.70000-01 -1.1G28M7) 
-2.50000-01 +1.8320H7i 
-2.00230-01 -2.00100-011 
-2.00230-01 +2.00100-011 
-1.33030-01 -7.20530-01i 
-1.33010-01 +7.20530-011 
-1. 50000-01 -1. 00000-01 i 
-1.50000-01 +1 .00000-01 i 
-1.40000-01 +1.40000-01) 
-1.40006-01 -1.40006-01) 
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'ON HOWS 



CLOSED-LOOP TIME RESPONSE 
MACH NO VS TIME FOR A RAMP MACH NO COMMAND 
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ALTITUDE 



TIME 


CLOSED-LOOP TIME RESPONSE 
ALTITUDE VS TIME FOR A RAMP MACH NO COMMAND 
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APPENDIX A- IV 

DESIRED EIGENVALUES AND EIGENVECTORS 
USED IN THE SYNTHESIS 


A-IV-1 


DESIRED EIGENVALUES: 


-2.00009-01 *2. 00009-01 i 
-2.00009-01 -2. 00009-01 i 
-1.50009-01 +1. 00009-01 i 
-1. 50009-01 -1. 00009-01 i 
-1.40009-01 +1. 40009-01 i 
-2.50009-01 *0.00009+00i 
-1.40009-01 -i.4000H>li 
-2.70009-01 *0.00009+00i 


DESIRED EIGENVECTORS: 


t 

p 

h 

h 

MEASUREMENTS 

a x M 

/M 

/ h 


c. 

0. 

0. 

0. 

99. 

99. 

99. 

0. 

V 

0. 

0. 

99. 

99. 

1. 

99. 

0. 

99. 

a 

0. 

0. 

99. 

99. 

0. 

0. 

0. 

99. 

q 

0. 

0. 

0. 

1. 

99. 

0. 

99. 

0. 

e 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

e 

99. 

1. 

0. 

99. 

0. 

0. 

0. 

0. 

p 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

r 

1. 

99. 

99. 

99. 

0. 

0. 

0. 

0. 


0. 

0. 

i. 

99. 

99. 

0. 

99. 

99. 

h 

0. 

0. 

0. 

0. 

99. 

1. 

0. 

A Thrust 

99. 

99. 

99. 

99. 
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APPENDIX B 


A DEMONSTRATION EXAMPLE FOR EXACT NONLINEAR 
CONTROLLER DESIGN WITH PRELINEARIZING TRANSFORMATION 
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j NjHomic r ion 

The theory of prel inoar 1 zi ng transf ormati ons for a class of nonlinear 
feedback system design has already achieved maturity. Pioneering work on 
these transformations for flight control problems was carried out at NASA 
Ames research center by G. Meyer and others [1-1A]. Though several papers 
are currently available, they deal with large dimension problems and conse- 
quently do not contain sufficient detail to serve as introductory example. 

The objective of this technical memo is to illustrate the application 
of these transformations with a simple but realistic example. The theory is 
excluded completely from the discussions since the transformation scheme for 
this example turns out to be direct. The performance of the exact nonlinear 
controller is compared with a gain scheduled linear perturbation controller 
to bring out the advantages of the former in clear detail. 

ILLUSTRATIVE EXAMPLE: 


Consider a two state variable model of a transport aircraft landing in 
the presence of winds as follows 


V 


nT-D 
mV cosY 


- ^tanY 


( 1 ) 


h - tanY 


( 2 ) 


here, V is airspeed, h altitude, T maximum thrust, D aerodynamic drag, 

Y flight-path angle, and n is the throttle setting. Y and n are the control 
variables in this model and down range is the independent variable. Now, a 
hypothetical thrust and drag characteristics are assumed as 


T - T q (V ♦ V w ) 

D - e~ Bh D (V ♦ V ) 2 
o w 


E-2 


and 


V - V (h) f V is the wind spend specified along flight path. 

WWW 

Note that the engine model is grossly incorrect at low speeds. 

It is required that the aircraft follow a given trajectory V R (t), 

h R ( 1 ) ' 

Since there are two state variables and two control variables in this 
problem, consider a transformation of the form 


X 1 “ V. x 2 “ h 


u 


1 


— -p— ■ — * $ tanY ; u 
mV cosY V 2 


tanY 


such that the transformed system is the form 


X 1 “ u 1 


X 2 “ U 2 


(3) 

(^) 


This model is in Brunovsky’s canonical form according to Meyer [10]. 
If the number controls are less than the number of states, a simple trans- 
formation such as the one given above is no longer possible. A transfor- 
mation is still feasible under certain weak conditions but will involve 
additional algebra. 

To consider the regulator problem first, a controller of the form 


r 1 

u i 


’ k 11 

r 

CM 

-X 


p— — < 

X 1 



T 

k 12 


— 

V 

U 2 


_ k 21 

k 22_^ 


- X 2. 


_ k 21 

k 22 J 


h 
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can !»' designed for the system (3), (M using any of the aviilaMo* techni- 
ques for linear system design. Note, however that the designed controller 
should satisfy physical constraints on u and 11 . 

Real controls can then be computed as 


tan 1 (k 21 V + k^h) 


n - [k,,V ♦ k, _h + •§ tanY] 


mV cosY D 


11 


12 


4 — 

T 


(5) 

( 6 ) 


Expressions (5) and (6) give the nonlinear regulator for the aircraft 
problem. A block diagram of the plant with regulator is given in Fig. 1. 


+ 


tan * (k^ j V + ^22^ 


[k, V + k ,h + ftany] 


trVcosy D 

T T 


Nonlinear 
Aircraft 
Model 

TF 

Disturbances 


Figure 1. Nonlinear Regulator with Aircraft 


The tracking problem is addressed next. Define the new state 
variables 


x - V-V R 
y - h-^ 




ORIGINAL PAGE IS 
OF POOR QUALITY 




X 


nT-P 

Ri Ji'OS) 


- ^ tarr>-C i 


y - tanY-C^ 


where 


h ■ C 
R 2 


using the transformations 


x 


1 


u 


u 


1 

2 


X, x 2 - y 

H T ‘ D - - i 

mV cosy V 

tanY-C 2 


tanY-C 1 


one has the following 


x 


2 


u 


u 


1 

2 


(7) 

( 8 ) 


which is in the Brunovsky's canonical form. 

As before, a feedback regulator can be synthesized for (7), (8) as 
shown below using linear design approaches. 


1 — — 

U 1 



k 12 


V 


\ 



m 





■ 



__ U 2_ 


f 

TT 

rv> 

k 22_ 


_ X 2 


^ k 21 

k„ 

22 


This controller can be transformed to original coordinates to give: 


Y - tan"* 1 l* 21 (V-V R ) 4 k 22 (h-h R ) 4 y 


(9) 


1 


E-5 



4 


( 10 ) 


n - fk n (V-V R ) * k 1 ? (h-h H ) « C, 


E 

v 


. , mV cosY 

tnnY] — 


D 

T 


Expressions (9) and (10) give the non] inear controller for tracking the 
commands V R (t) and h R (t). If a tighter control is desired, additional inte- 
gral feedbacks can be incorporated with very little difficulty. 


NONLINEAR CONTROLLER EVALUATION 


To evaluate the nonlinear controller (9) and (10), a feedback gain 
matrix was computed first using linear quadratic regulator theory with the 
following state and control weighting matrices. 


R 

0 


R 

" 5000 

0 

XX 0 

100 

t 

UU 

0 

10® 


These computations were carried out with the hypothetical data 
T D 

— - 0.0*1025, — - 0.0136238 x 10 . 

m m 

To serve as a standard for comparison, a gain scheduled linear pertur- 
bation controller was designed with the same data. Four design conditions 
were used which are given in Table 1. 


TABLE 1 . LINEARIZATION POINTS FOR THE ILLUSTRATIVE 
EXAMPLE 


Range 

feet 

h R* 

h 

feet 

V 

V 

feet/s 

C 2’ 

• 

h 

V 

• 

V 

0 

10000 

350 

“0.0*1667 

“0.000*167 

150000 

3000 

280 

“0.06 

“0.0016 

200000 

0.0 

200 

0.0 

“0.02 

205000 

0.0 

60 

0.0 

0.0 
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Note tint, the linearization points were on the commanded trajectory. 

Two simulations are next setup. The first one implementing the non- 
linear controller (9) and (10) and the second one using the gain scheduled 
linear perturbation controller. A disturbance in the form a wind shear was 
introduced in the simulations for additional realism, (fig. 2). 

Figures 2-4 give the results. Note that the nonlinear controller has 
a slightly superior tracking performance. The real advantage of using this 
controller, however, is that it performs well using one set of gains 
throughout the operating region while the linear perturbation controller 
requires a gain schedule to achieve acomparable performance. 

CONCLUSIONS 

A simple illustrative example was discussed in this technical memo to 
demonstrate the power of the prelinearization transformations in 
synthesizing exact nonlinear controllers. Comparisons were made against a 
gain scheduled linear perturbation controllers to show relative advantages. 


12000 


10000 


BDDO 


1 2 


Figure 2 


34567 & StO 

V»INDSPEED(FT/S) 


Wind Shear Used in Simulations 
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Figure 3. Comparison of the Performance of Linear and Nonlinear Con- 
trollers: Altitude vs. Range 


ALTITUDE 


i 



RANGE 


Figure *J. Comparison of the Performance of Linear and Nonlinear 

Controllers in the Vicinity of Touchdown: Altitude vs. Range 

A: Nonlinear Controller 

B: Linear Perturbation Controller 

C: Commanded Trajectory 
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AIRSPEED 



Figure 5. Comparison of the Performance of Linear and Nonlinear Con- 
trollers: Airspeed vs. Range 

A: Nonlinear Controller 

B: Commanded Trajectory 

C: Linear Perturbation Controller 
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APPENDIX C 


FORTRAN PROGRAM FOR MANEUVER MODELING 
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3-D INTERPOLATION : RIGHT HANDED TRIAD DATA 1/1/ 

X-ALTITUDE C 
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SUBROUTI NE ATMO ( ALTI TU , VSOUND , DENSTY , VI SCO , DVSDH , DDDH , DVDH) 
IMPLICIT REAL'S (A-H,0-Z) 

DIMENSION ALT (9) , SPEED (9) ,DENSIT(9) ,VIS(9) 


APPENDIX D 


LINEAR TIME VARYING SIMULATION 
IN SYSTEM BUILD 
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LINEAR TIME VARYING SIMULATION IN SYSTEM BUILD 


Linear time varying simulation can be built around a Fortran block 
evaluating the matrix equation 

y = A(t)X + B(t)U (D. 1 ) 

where y e I*’, X c B n , U e E m . A and B are time varying matrices of com- 
patible dimensions. To use this Fortran block, the A and B matrix 
elements must be arranged in the following tabular form and transferred to 

Matrix stack, 
x 



t 

0 

A 11 

A 21 

A , 
nl 

A 12 

A A 

22' ' ' n2' * ’ 

A 

nn 

B n 

B 21 

B n1 

... B 

nm 



s 

A 1 1 

A 21 •*' 

A , 
nl 

A 1 2 

A 22’ " A n2 

A 

nn 

B n 

B 21 

B n1 

... B 

nm 


Table = 

l 2 

A 1 1 
• 

A 

21 

• 

A i 

nl 

A 12 

A 22* -A n2 

• • 

A 

nn 

• 

B n 

B 21 

• 

B n1 

• 

... B 

nm 

• 

(D.2 


W * 

• 

% 

♦ 

. 

• % 

• 

• 

• 

• 

* ^ 



independent 

variable 

To build in the Fortran block, two distinct steps are necessary. The first 
one consists of building a dummy superblock with the name TABLE(N) , where N 
could be any number from 0 to 9. The Table D.2 is built as a linear 
interpolation table block in this superblock, say at relative location M. 

Next, the superblock(s ) including the time varying Fortran block may 
be built with the following responses for the queries 
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Dimension of RP & IP: 0,2 (no real parameters, two real parameters) 

enter IP : N,M (Dummy superblock Table number, relative 

location) 

Note that to complete the time varying simulation, one has to define 
an independent string of l integrators. A schematic flow chart for building 
the time varying linear system with this FORTRAN block is given in Figure 

D.l. 

TM 

The SYSTEM_BUILD block diagram for linear time varying simulation 
for the aircraft example discussed in this Appendix is given in Figure D.l. 
The first block in this diagram computes the derivatives of the states using 
the Fortran block. These are then integrated in block 2. Block 3 
essentially routes various inputs and outputs. The reference trajectory is 
generated in Block H while the wind is generated in Block 5. Block 6 imple- 
ments the time varying gains or gain schedules, again using the Fortran 
block . 
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N: Any number from 
0 to 9 


M: Any number from 
t to 6 


Figure D.l. Flow Chart for Building Linear Time Varying Simula ’or. on 
SYSTEM BUILD 
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MINIMUM ERROR EXCITATION OUTPUT FEEDBACK DESIGNS 


This appendix gives the performance index specification for the 
control laws designed at each flight condition. Together with the algorithm 
description of section *1.2.2, these performance indeces completely specify 
how the 30 sets of output gains can be reproduced. 


LQR Performance Index 


J = J 0 e ^ at (x^HA Hx + u^R u)dt, with 

JO y uu 

. r 1 1 1 1 1 1 

y gL (50) 2 * (0.01) 2 ’ (0.001) 2 ’ (0.001) 2 ’ (0.001) 2 ’ ( 0 . 001 ) 2 * 

111 11 
(0.001) 2 ’ (0.001) 2 ’ (0.001) 2 ’ ( 0 . 005) 2 ’ (0. 005) 2 ’ 

3 ] 

(0. 005 ) 2 


State weighting matrix, R^ x 


k H T A H, 

y 


here: k is a scalar providing relative state-control weights 

H is the system output matrix 


The control weighting matrix R u|j = diag 2 (B) 

Table D-1 lists the scalar weight k and the state weighting vector B as a 
function of the design flight condition. 
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TABLE E-1. PERFORMANCE INDEX SPECIFICATION FOR EACH FLIfflT CONDITION 
(Subscript 1 corresponds to the design with all controls, 
subscript 2 corresponds to the design without throttle 
control . ) 


ALTITUDE AND MACH 

LOAD h - 10k h « 20k h - 30k h - HOk h - 50k 

FACTOR M - 0.8 M - 1 .0 M-1.2 M - 1 . M M-1.8 


B,- [0.1 111] B] 
k, - 0.002 k, 

B 2 - tl 1 1] B 2 

k 2 - 0.002 k 2 


B,- [5111] B, 

k, - 0.001 k, 

B 2 - [0.1 0.1 1] B 2 

k 2 » 0.0005 k 2 


B,« [5111] B, 

k, - 0.001 k ! 

B 2 - [0.1 0.1 1] B 2 

k 2 - 0.002 k 2 


[ 0.1 1 1 1 ] 
0.002 


B, - [0.1 1 1 1] 
k, - 0.001 


B, - [0.1 1 1 1] 
k, - 0.001 


B , = [0.1 1 1 1 ] 
k, - 0.001 


[1 1 1] 
0.002 


[5111] 

0.001 


B 2 - [1 1 1] 
k 2 - 0.005 


B, - [5 1 1 1] 
k, - 0.001 


B 2 

k 2 


B, 

k 1 


-[111] 

- 0.01 


-[5111] 

- 0.001 


B 

k 


B 

k 


2 - [1 1 1 ] 

2 - 0.01 


i - [5 1 1 1] 

, - 0.002 


[0.1 0.1 1] B 2 - [0.1 0.1 1] B 2 
0.001 k 2 - 0.001 k 2 


[5111] B, - [5 1 1 1] B, 

0.001 k, - 0.0005 k, 


[ 0.1 0.1 1 ] 
0.01 


[51 1 1 ] 
0.0005 


B 

k 


B 

k 


2 - [ 0.1 0.1 1 ] 

2 - 0.00001 


,-[ 5111 ] 

, - 0.00012 


[0.1 0.1 1] B 2 - [0.1 0.1 1] B 2 
0.002 k 2 - 0.002 k 2 


[0.1 0.1 1] b 2 

0.002 k 2 


[ 0.1 0.1 1 ] 
0.00001 


Output Feedback Gains 


The output feedback Controller has the form 


u = Cy 


where 


u - [ST se ap sa ap sr ap : 


y = [6h 6M 6a 5 Y 641 6B 6p 6q fir /6a / 6<f 


Flight 

condition: 

: h-lOk 

ft, M=0.8 

, n=1 



(a) 

Throttle free 






COLUMNS 1 

6.1491D-04 
-2.1757D-02 
“8 . 5700D-05 
*3.42160-06 

THRU 6 
1 .9996D+03 
9.3094D+00 
2.7476D-02 
1.3714D-03 

1 . 2553D+01 
-4.81150+01 
-2 .50240-01 
-9.5302D-03 

3.91280+03 
-3.81090+01 
-1 . 2745D-01 
-5.2506D-03 

3. 0025D-02 
-5.8576D-01 
5.8288D+01 
6 . 184 ID- 01 

7.01120-02 

-3.0951D-:i 

-9. 71120+01 
3. 3190D+00 


COLUMNS 7 
-3 . 5007D-02 
-3.7222D-01 • 

2 . 3275D+01 • 

8.1391D-01 ■ 

THRU 12 
2.78330+00 
-2 . 5667D+01 
-1.5444D-01 
-7.54310-03 

-6.43620-03 
1.61860-01 
2 . 3048D+01 
-4.39780+00 

-5.88240+03 

-6.4055D+01 

-2.4475D-01 

-1.0375D-02 

1.8056D-02 

-3.17820-01 

3.2259D+01 

3.7226D-01 

3. 61500-03 
-5. 41630-02 
5.48720+00 
6. 35600-02 

(b) 

Fixed throttle 






COLUMNS 

1 THRU 6 






0. OOOOD+OO 
-6 . 9831D-01 
-1.9239D-03 
-1.57750-05 

0 . OOOOD+OO 
-1 . 7886D+04 
-4 . 9332D+01 
-4.0474D-01 

0. 0000D+00 
-9 .43240+01 
-2.6985D-01 
-2 . 2771D-03 

0. 0000D+00 
-2 . 5680D+02 
-6 . 9729D-01 
-5.6596D-03 

0. OOOOD+OO 
-1 . 9276D-01 
5.8289D+01 
6. 1835D-01 

0. OOOOD+Oo 
-2.7120D-01 
-9 . 7113D+01 
3 . 3189D+00 


COLUMNS 7 
0 . OOOOD+OO 
-1.5857D-01 
2 . 3276D+01 
8 . 1388D-01 

THRU 12 
0. 0000D+00 
-2.07620+01 
-1 . 2968D-01 
-9 .42640-04 

0. 0000D+00 
2 . 0846D-01 
2 . 3048D+01 
-4 . 3978D+00 

0. 0000D+00 
0. OOOOD+OO 
0. OOOOD+OO 
0. OOOOD+OO 

0. OOOOD+OO 
-9.4988D-02 
3 . 2260D+01 
3.7223D-01 

0. OOOOD+OO 
-1.5752D-02 
5.4873D+00 
6 . 3555D-02 
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Flight Condition: h=1 Ok ft, M-0.8, n=2 


(a) Throttle free 



COLUMNS 
2.9627D-02 
-B. 29390-02 
-6.9160D-03 
-3.5156D-04 

1 THRU 6 
8.1168D+02 
-1 .71090+03 
-3.03150*02 
-6.7290D+01 

-2.85920-02 
-4.34570+01 
5.38340+00 
3. 0272D+00 

-7.6192D-01 
-4 . 2727D+01 
3.8716D+01 
1.6798D+01 

-1 . 3855D-Q1 
-3 . 0855D+00 
4. 1267D+01 
-9.9O49D-01 

-3.80650-01 
-5 . 8524D-01 
-4.6182D+01 
1 . 2926D+01 


COLUMNS 
-8.1273D-03 
-1 .62310+00 
1.7513D+01 
-1.5747D-01 

7 THRU 12 
-3 . 1193D-01 
-1.0912D+01 
1.8824D+00 
7.68610-02 

6.5570D-02 
-2 . 1829D+00 
-3 . 1690D+00 
-3. 2354D+00 

7.7035D+01 
-1.9868D+02 
-4. 67120+01 
-1.4266D+01 

-2 . 8132D-01 
-6 .19470-01 
2.49010+01 
7 .70210-01 

-8.7290D-02 
-5. 45150-01 
4.06800+00 
-7.0858D-03 

(b) 

Fixed throttle 






COLUMNS 
0 . OOOOD+OO 
-9.4567D-01 
- 5. 452 0D- 02 
-1.8692D-02 

1 THRU 6 
0 . OOOOD+OO 
-2.25540+04 
-1.50960+03 
-5. 2788D+02 

O.OOOOD+OO 
-1. 75410+02 
3.5769D+00 
2-5422D+00 

O.OOOOD+OO 
-3.77940+02 
2 . 0487D+01 
5.5214D+00 

O.OOOOD+OO 
8 . 7822D+01 
1 . 7270D+02 
6.65570+00 

o.oooc: -dc 
-1.38090+02 
-6.21750+01 
5.6499D+00 


COLUMNS 
O.OOOOD+OO 
4 . 09910+01 
8.3971D+01 
3 . 9233D+00 

7 THRU 12 
0 . OOOOD+OO 
-1.26210+02 
1. 74010*00 
-3.6830D+00 

O.OOOOD+OO 
9. 67030+01 
-1.99410+00 
2.7026D+00 

O.OOOOD+OO 
-2 . 1918D+03 
-1.523BD+02 
-5.4415D+01 

O.OOOOD+OO 

3.4846D+01 

9.8190D+01 

4.51160+00 

O.OOOOD+OO 
6. 09210-00 
1.6852D+01 
7 . 3470D-01 

Flight Condition: 

h=1 Ok ft 

, M-0.8, n 




(a) 

Throttle free 






COLUMNS 
2 . 27940-03 
-4.3918D-02 
-3. 3214D-04 
2.01980-03 

1 THRU 6 
7.1210D+01 
-9.48D6D+02 
-2 .5846D+02 
-6.1265D*01 

4.02410-01 
-4.5341D*01 
2 . 8570D+00 
1 . 0238D+00 

-1. 15100+00 
-8.34750+00 
5.7110D+01 
2. 22730+01 

-6 . 1276D-02 
-1.6230D+01 
3.9784D+01 
-2 . 9065D+00 

-2 .43840-01 
1.4902D+01 
-4.6713D+01 
1 . 5259D+01 


COLUMNS 
9.7857D-03 
-1.7379D-01 
1.6825D+01 
-6 . 3254D-01 

7 THRU 12 
-8.7439D-02 
-7.5154D+0C 
2 . 34310+00 
2.44540-01 

1.9166D-02 
-1 .85270+00 
-3 . 2140D+00 
-3 . 4999D+00 

1.8054D+01 
-2 . 8085D+02 
-7 . 6799D+01 
-2.00670+01 

-3 . 3040D- 01 
7.6036D+00 
2 . 9976D+01 
2 . 14670+00 

-1 . 2460D-01 
6 . 2348D- 02 
5. 17820+ CC 
2.4978D-01 

(b) 

Fixed throttle 






C "LUMNS 
0.0D00D+00 
-2.4868D-01 
-1.3103D-02 
1.9612D-03 

1 THRU 6 
O.OOOOD+OO 
-4. 54800+03 
-2 . 0613D+02 
-3.6559D+01 

O.OOOOD+OO 
-1 . 9039D+02 
1 . 2893D+00 
1.6883D-01 

O.OOOOD+OO 

-1.8448D+02 

1.7878D+01 

1.4573D+01 

O.OOOOD+OO 
5 . 2790D+01 
2 .63140+02 
2.4480D+00 

0 .O0CCD«?0 
-3.5337D- 01 
-6.6656D+01 
8. 1336D+00 


COLUMNS 
0 . 0000D+00 
2 .76150+01 
1 . 3028D+02 
8 . 04870-01 

7 THRU 12 
O.OOOOD+OO 
- 1 . 6666D+ 0 2 
4.3779D+00 
-6 .89810-01 

O.OOOOD+OO 
4. 0532D+01 
-2. 27010+00 

-2.3050-;" 

O.OOOOD+OO 

-1.3112D+03 

-5.59570+01 

-1.23460-01 

O.OOOOD+OO 
2.18150+01 
1.51310+02 
3 . 0843D+00 

O.OOOOD+OO 

3.1686D+00 

2.6108D+01 

4.93320-01 


E-5 



Flight Condition: h-20k ft, M=1.0, n»1 


(a) Throttle free 



COLUMNS 

-3 . 0028D-02 
-2.4707D-02 
“6 . 2268D-05 
-1.5270D-06 

1 THRU 6 

1 .4232D+03 
3 . 9892D+00 
1.0190D-02 
3 . 3189D-04 

-1.4150D+02 
-I . 3192D+02 
-2.5252D-01 
-6 . 0395D-03 

1 . 2419D*04 
2 . 1845D+01 
6.6599D-02 
1.8372D-03 

1.5656D-01 
-2.2288D-01 
5.6690D+01 
1 . 0355D*00 

2.5076D-01 
-3 . 9225D-01 
-1 . 0369D*02 
3 .7093D+00 


COLUMNS 

7 THRU 12 






*3 . 1757D-02 
-1.9513D-01 
2 . 29840*01 
1.0205D+00 

-9.4699D+00 
-2.38180+01 
-1 . 3557D-01 
-4. 098 ID- 03 

“2.8527D-02 
2 . 9430D-01 
2 .5598D+01 
-4.4825D+00 

-1.7748D*04 

-1.8270D*02 

-4.6193D-01 

-1.2099D-02 

8.6553D-02 
-1 . 2200D-01 
3. 1368D- Cl 
6.07500-01 

1.6153D-02 
-2.09720-02 
5 . 3333D-00 
1 . 0334D-01 

(b) 

Fixed throttle 






COLUMNS 
0 . OOOOD+OO 
-1.4623D-01 
-3.4606D-04 
-8 . 5753D-06 

1 THRU 6 
0 . 0000D+00 
-3 . 0862D+03 
-7 . 3324D+00 
-1.8118D-01 

0 .0000D+00 
-1 . 3772D+02 
-3 . 2365D-01 
-7.6613D-03 

0 . 00000*00 
-2 . 5546D*02 
-5.9666D-01 
-1.4886D-02 

0 . 00O0D*00 
-2 . 3182D-01 
5.6691D*01 
1 . 0355D*00 

0 .OOOOD+OO 
-3.56610-01 
-1 . 0369D*02 
3.70930-00 


COLUMNS 
0 .OOOOD+OO 
-1.89980-01 
2 . 2985D+01 
1 . 0205D+00 

7 THRU 12 
0 . OOOOD+OO 
-2 . 11460+01 
-1.4990D-01 
-4 . 2346D-03 

0 . 0000D+00 
2.7138D-01 
2.5599D+01 
-4.4824D+00 

0 . 00000*00 
0.00000*00 
0 . OOOOD* jO 
O.OOOOD*CO 

0 . 00000*00 
-1. 1843D-01 
3 . 1368P~C1 
6.074^-Cl 

0 . OOOOD-C 0 
-2. 00580- '2 
5.33341- 
1.03341 

Flight Condition: h«20k 

ft, M* 1 .0 

, n-2 



(a) 

T nrottle 

free 






COLUMNS 
6 . 9042D-03 
-8 . 0054D-02 
-3.1811D-03 
7 . 3580D-04 

1 THRU 6 
1.9549D+02 
-1.6813D+03 
-I . 9103D+02 
-2 . 9358D+01 

1 .5364D+00 
-9 .5674D+01 
7 .1812D+00 
3.8067D+00 

-8 .59170-01 
-3.7453D*01 
4 . 0158D*01 
1.8625D*01 

1 . 1123D-01 
-1.7710D*01 
4 . 1578D+01 
-1 . 00290*00 

-4.6072D-01 
1. 1887D+01 
-4.8085D+01 
1 . 2088D+01 


COLUMNS 
1 . 1787D-01 
~1 . 200 3D* 01 
1 .7896D+01 
8 . 7576D-03 

7 THRU 12 
-3 . 2178D-02 
“1. 7947D+01 
1 . 8398D+00 
-2.7686D-02 

1.3241D-01 
-8.07140*00 
-2 .7051D+00 
*3 . 0579D+00 

1.8588D+01 
-2 . 1281D+02 
-3.47230*01 
-1.08970*01 

-1.5842D-01 
-7.5076D+00 
2.4838D*01 
9 . 1879D-01 

-5.7223D-02 
-1.4183D+00 
4 . 0935D+00 
-1.1976D-02 

(b) 

Fixed throttle 






COLUMNS 

1 THRU 6 






0. OOOOD+OO 
*2 .5821D-01 
-3. 338 ID- 03 
-1.4984D-03 

0.00000+00 
-5 . 1355D+03 
-1 .8932D+02 
-9 . 3411D+01 

0 . 00000*00 
-1.8856D+02 
8.86310*00 
1 . 1209D*00 

0 . 0000D*00 
-2 .4846D+02 
3 . 0696D+01 
1. 3272D*01 

0 . 00000*00 
2 . 3176D+01 
2 . 0976D+02 
2.82450*00 

0. OOOOD+OO 
-5.6307D+01 
-5.8322D+01 
7. 2034D+00 


COLUMNS 
0 . 0000D+00 
9 .70170+00 
1 . 0345D+02 
1 . 2950D+00 

7 THRU 12 
0 . OOOOD+OO 
“1 . 1540D+02 
7.8005D+00 
“1 . 1666D+00 

0 . OOOOD+OO 
2.1620D*01 
-3.3573D+00 
“2.2116D*00 

0 . 0000D*00 
-5.6373D+02 
-2.34740*01 
-1.4181D*01 

0 . 00000*00 
6.80680*00 
1.19370*02 
2.49880*00 

0. OOOOD+OO 
1 . 0118D+00 
2 . 0505D+01 
3.6095D-01 


E-6 


Flight Condition: h = 20k ft, M=1.0, n=^ 


(a) Throttle free 



COLUMNS 
-2 - 9947D- 05 
-2.4350D-02 
8 . 9437D- 03 
4 . 1491D-03 

1 THRU 6 
3.6599D+00 
-1.6104D+02 
-5 . 3990D+01 
-1 . 8467D+01 

5 . 5165D-01 
-9 . 3784D+01 
5 . 8315D+00 
5 . 9942D-01 

-4 . 3026D-01 
-5 . 3770D+01 
6 . 7524D+01 
2.9591D+01 

1 . 0464D-01 
9 . 3674D+ 0 0 
3.8510D+01 
-3.6880D+00 

-1.4545D-01 
-2. 07840+01 
-4.8181D+01 
1.6136D+01 


COLUMNS 
1.7035D-02 
3 .7699D-01 
1 .6541D+01 
-1 . 970 3D-*- 00 

7 THRU 12 
1 . 0704D- 01 
-2 .8745D+01 
1 .7075D+00 
-3.2254D-01 

1 . 9226D-02 
1 . 1755D+0Q 
-3 .4304D+00 
-4 . 4090D+00 

1 .4165D+00 
-9.7631D+01 
-2 . 9706D+01 
-1 . 1588D+01 

-3.8941D-02 
4.2783D-01 
2 . 9120D+01 
9 . 5902D-01 

-1.7625D-02 
2.9496D-01 
5 . 0757D+00 
1 . 5582D-01 

(b) 

Fixed throttle 






COLUMNS 
0 . 0Q00D+00 
-6 . 01330-02 
1.3187D-03 
3.5796D-03 

1 THRU 6 

0. 0000D+00 
-4 . 3196D+02 
9 . 5875D+00 
-4.9691D-01 

0 . OOOOD+OO 
-1.5452D+02 
1 . 5789D+01 
1 . 9333D+00 

0. OOOOD+OO 
-1 . 3945D+02 
3 . 1355D+01 
1 . 9762D+01 

0. 0000D+00 
1 . 1041D+01 
2.5958D+02 
3.1445D-01 

0.0000D+0C 
-4.5528D+01 
-6.5193D+01 
7 . 3005D+00 


COLUMNS 
0 . OOOOD+OO 
3 . 4082D+00 
1.2898D+02 
-4 7916D-01 

7 THRU 12 
0. OOOOD+OO 
-1 . 4477D+02 
1 . 7756D+01 
2.1208D+OC 

0 . OOOOD+OO 
3.1317D+00 
-3 . 461 0D+ 00 
-4 . 14330*00 

0 , 0000D+0Q 
-2 . 7722D+02 
9 . 5732D-00 
-1. 0696D-00 

0. 0000D+00 
-5 . 3801D+00 
1 . 4989D+02 
1 . 6873D+00 

0. OOOOD+OO 
-1.25550+00 
2 . 5871D+01 
2. 90840-01 

Flight 

condition: h - 30 k ft, M- 1 

.2, n® 1 



v a ) 

Throttle : 

free 






COLUMNS 

-4.6Q88D-02 

-5.7504D-02 

-4.0537D-05 

-7.3106D-08 

1 THRU 6 
1.6086D+03 
-2 .44810+01 
1.7568D-02 
1.2564D-04 

-4 . 0862D+01 
-1 . 9418D+02 
-2 . 066 2D- 01 
-2 . 9104D-04 

7 . 7420D+03 
5 . 1413D+01 
1.46920-01 
2.96620-04 

5 . 0626D-01 
1 . 9622D+00 
4 . 2962D+01 
-2. 29940-01 

2 . 2822D-01 
4.4702D-01 
-7. 5834D+01 
8.08500+00 


COLUMNS 
2 . 5531D- 01 
1 . 0298D+00 
1 .67290+01 
6 .75630-01 

7 THRU 12 
1 . 96770+01 
-5 . 9481D+00 
-7.8761D-02 
-1 . 0035D- 03 

-8.0973D-03 
1.66250-01 
1 . 1728D+01 

-5.5422D+C0 

-6 . 7504D+03 
-2.3518D+02 
-3. 0407D-01 
-7.55420-04 

2 . 9316D-01 
1.1180D+00 
2 . 3583D+01 

-4.14570-02 

5.1002D-02 
1 . 9185D-01 
4 . 0053D+00 
-4.9060D-03 

(b) 

Fixed Throttle 






COLUMNS 
0 .OOOOD+OO 
-2 . 2035D-01 
-3.9297D-04 
-9.48180-06 

1 THRU 6 
0.000 0D+ 0 0 
-4 . 5955D+03 
-8 . 2035D+00 
-1.9769D-01 

0 . OOOOD+OO 
-2 .4374D+02 
-3 . 9043D-01 
-9 . 3204D-03 

0. OOOOD+OO 
-4.8065D+02 
-8.4043D-01 
-2 . 0343D-02 

0. OOOOD+OO 
1 .45410+00 
7.7573D+01 
1.90120+00 

0. OOOOD+OO 
-1.1746D-01 
-9 . 1042D+01 
2.6540D+00 


COLUMNS 
0 .OOOOD+OO 
7 . 5137D- 01 
3 . 4197D+01 
1. 3758D+00 

7 THRU 12 
0. OOOOD+OO 
-1.7241D+01 
-1. 58160-01 
-4.1546D-03 

0 . OOOOD+OO 
2 . 9844D-01 
2.8032D+01 
-5. 14480-00 

0 . 0000D+00 
0 . 0000D+00 
0 . OOOOD+OO 
0 . OOOOD+OO 

0. OOOOD+OO 
8.2581D-01 
4 . 3223D+01 
1 . 0546D+00 

0. OOOOD+OO 
1 .4191D-01 
7 . 3793D+00 
1 . 7917D-01 


E-7 



Flight 

Condition: 

h=30k ft 

, K-1.2, 

n=2 



(a) 

Throttle free 






COLUMNS 
6 . 6424D- 03 
-7.9027D-02 
-2 . 569QD- 03 
1.9215D-03 

1 THRU 6 

2 . 0838D+02 
-1.9271D+03 
-2 . 0704D+02 
-2 . 3017D+00 

2 . 4550D+00 
-1 . 7763D+02 
9 . 0405D+00 
6.6459D+00 

-4.6783D+00 
3 . 9364D+01 
4.77580+01 
1.90070+01 

2 .77070+00 
-1 . 1536D+02 
4.4819D+01 
2. 16340+00 

-2 . 5482D+0C 
5.3945D+0I 
-3.83430+01 
7. 74380+ OC 


COLUMNS 

7 THRU 12 






1 . 3539D+00 
-8 . 0765D+01 
1.8876D+01 
2.45010+00 

1 . 8265D-02 
-1 . 9153D+01 
1. 34490+00 
5.2374D-01 

1.0172D+00 
-3.15580+01 
-3 . 4281D+00 
-1 . 5152D+00 

4. 1165D+01 
-4 . 4553D+02 
-5 . 0063D+01 
-3 , 9816D+00 

9. 1485D-01 
-5 . 0057D+01 
2.6362D+01 
3 . 1590D+00 

1.2343D-01 
-7 . 8047D+00 
4. 3945D+00 
2.5665D-01 

(b) 

Fixed throttle 






COLUMNS 
0 . 0000D+00 
-2 . 7603D-01 
-1 . 0980D- 02 
-2 . 3Q01D- 03 

1 THRU 6 
0 .OOOOD+OO 
-5 . 6529D+ 03 
-4. 1304D+02 
-1 . 2616D+02 

0. OOOOD+OO 
-2 . 9060D+02 
3.7056D+00 
1.4133D+00 

0. OOOOD+OO 
-3.6612D+02 
4.25480+01 
1.8791D+01 

0. OOOOD+OO 
5 . 2588D+01 
2 . 1125D+02 
3 . 0150D+00 

0. OOOOD+OO 
-8.87930+01 
-3 . 8144D+01 
9 . 2470D+0 0 


COLUMNS 
0 . 0000D+00 
1 . 8824D+ 01 
1 . 0450D+02 
2.2519D+00 

7 THRU 12 
0 .0000D+00 
-9. 40440+01 
4. 14140+00 
-2.2669D+0C 

0. OOOOD+OO 
4 . 0305D+01 
-7 . 5652D+00 
-2.C4480+CC 

0. OOOOD+OO 
-1.1392D+03 
-8.79910+01 
-2.92380+01 

0. OOOOD+OO 
2 . 0379D+01 
1 , 2094D+02 
3. 50120+00 

0 .OOOOD+OO 
3 . 5850D+00 
2 . 0762D+01 
4 . 7414D-01 

Flight 

Conditions 

: h=30k ft, M=1 .2, 

n=*J 



(a) 

Throttle free 






COLUMNS 
4. 1769D-04 
-6.6101D-02 
6 . 4192D-03 
4 . 5633D- 03 

1 THRU 6 
3.6498D+01 
-1 . 0172D+ 03 
-2 .86740+02 
-4 . 2698D+01 

1.4851D+00 
-2 . 1118D+02 
4 .56540+00 
4 . 9646D+00 

-3.4841D+00 
-1 . 4182D+02 
1 . 1561D+02 
4 . 3199D+01 

2. 10420+00 
-1 . 1462D+01 
2.7896D+01 
-6.6763D+00 

-1 . 3394D+0C 
-7.41840+01 
-2.6873D+01 
2 .01810+01 


COLUMNS 
6.36310-01 
-3 . 1037D+ 01 
1 . 0313D+01 
-2.4702D+00 

7 THRU 12 
4.8008D-02 
-2 . 6264D+01 
8.3152D-01 
3 . 2870D-01 

3. 3416D-01 
1.77620+00 
-4 . 9400D+00 
-4.1412D+00 

1.4069D+01 
-4.07610+02 
-1 . 1891D+02 
-2 . 0022D+01 

5.33120-01 
3 . 2469D+00 
2.71990+01 
9. 2925D-01 

9 . 1566D-02 
7 . 1053D+0 0 
4 . 8492D+00 
-2 . 3495D- 02 

(b) 

Fixed throttle 






COLUMNS 

1 THRU 6 






0 . 0000D+00 
-1.4369D-01 
-1.8925D-03 
2. 3785D-03 

0. OOOOD+OO 
-2.6046D+03 
-1 . 94040+02 
-4.3929D+01 

0.00000+00 
-2 . 7650D+02 
2.7848D+00 
-1 . 3042D+00 

0. OOOOD+OO 
-2.2663D+0 2 
6.79970+01 
2.6549D+01 

0. OOOOD+OO 
1.87780+01 
2 . 8567D+02 
-4. 3513D+00 

0.00000+00 
-6.04320+01 
-3. 9003D+01 
1. 0744D+01 


COLUMNS 
0.00000+00 
3 . 3267D+00 
1. 4358D+02 
-2.1208D+00 

7 THRU 12 
0. OOOOD+OO 
-1 . 3238D+02 
5.40950+00 
-8.7326D-01 

0. OOOOD+OO 
2.55970+01 
-1 . 5841D+01 
-4 . 3539D+00 

0. OOOOD+OO 
-1 . 0794D+03 
-7 . 9803D+01 
-1 . 9425D+01 

0. OOOOD+OO 
-2.48470+00 
1.6880D+02 
2 .88940-01 

0 . OOOOD+OC 
-7.81230-01 
2 . 9097D+01 
4. 05360-0' 


E-8 


Flight Condition: h= i )0k ft, m=1.*), n=1 

(a) Throttle free 


COLUMNS 1 THRU 6 



-4.0290D-02 
-3 . 3645D-02 
-2.2849D-05 
7 . 4382D-07 

1 . 4738D+03 
-3 .4998D+01 
1 . 0441D- 02 
-2.6323D-04 

-7 . 8132D+01 
-2 . 2104D+02 
-1.9676D-01 
6 . 3096D-03 

1 . 2733D+04 
2 . 6056D+02 
4 . 0965D- 01 
-1 .31870-02 

2 . 5962D+00 
8.6121D+00 
4 . 4062D+01 
-1 . 5759D+00 

4.70390-01 
1.64250+00 
-5 . 9481D+01 
1 . 1276D+01 


COLUMNS 
1.5059D+00 
4 . 9056D+00 
1 . 8275D+01 
1 . 2100D-01 

7 THRU 12 
2 . 1803D+01 
-1 . 1218D+00 
-7 . 3752D-02 
1.6192D-03 

1 . 5995D-01 
4 . 6368D- 01 
6 . 8008D+00 
-7 . 5773D+00 

-7.5485D+03 
-3 . 34570+02 
-3.9080D-01 
1 . 2509D-02 

1.4565D+00 

4.7969D+00 

2.4078D+01 

-7.5866D-01 

2.4843D-01 

6.15300-01 

4.08750+00 

-1.2434D-01 

(b) 

Fixed throttle 






COLUMNS 
0 . 00000+00 
“1 . 7781D-01 
-3.6355D-04 
-5.8903D-06 

1 THRU 6 
0 . 0000D+00 
“1 . 0549D+04 
*2 . 1593D+01 
-3.4989D-01 

0 .00000+00 
-1.96240+02 
-3 .6862D-01 
-5 . 9813D-03 

0. OOOOD+OO 
-9. 2410D+02 
-1.87570+00 
-3 . 0397D- 02 

0.00000+00 
1.11140+01 
1 . 0197D+02 
2.25520+00 

0. OOOOD+OO 
6.75300-01 
-7.48400+01 
6 . 3920D+00 


COLUMNS 
0. OOOOD+OO 
6 . 0610D+00 
4.77440*01 
1 . 7109D+ 00 

7 THRU 12 
0 . OOOOD+OO 
-1 .7085D+01 
-1 .60070-01 
-2.67490-03 

0. OOOOD+OO 
9 . 1804D-02 
2 . 1718D+01 
-8 . 6307D+00 

0. OOOOD+OO 
0 .OOOOD+OO 
0. OOOOD+OO 
0 . OOOOD+OO 

0. OOOOD+OO 
6 . 2255D+00 
5 . 6830D+01 
1 . 2089D+ 00 

0. OOOOD+OO 
1 . 0664D+00 
9 . 7189D+00 
2 . 0712D-01 

Fli 

ght Condition: h-40k ft, 

IV 2 



(a) 

Throttle free 






COLUMNS 
6.8518D-03 
-4 .45510-02 
-9.5491D-04 
2.1000D-03 

1 THRU 6 
5 , 06830+02 
-2 . 32510+03 
-3 . 2593D+02 
-1 .5805D+01 

4 . 0818D+00 
-1 . 8765D+02 
9.99910+00 
9 . 0181D+00 

-5.6239D+00 
-7.07420+01 
7 . 9025D+ 01 
4.9655D+01 

3 . 0230D+00 
-5.8678D+01 
4. 3050D+01 
-4.8933D+00 

-1.7841D+00 
1.4762D+00 
-3 . 0872D+01 
1.5230D+01 


COLUMNS . 
1.4281D+00 
-6 . 50570*01 
1. 80560*01 
-6.8358D-01 

7 THRU 12 
-1 . 3268D-01 
-1 . 5583D+01 
1 . 2552D+00 
5.66240-01 

7 .72010-01 
-1.2903D+01 
-4 . 8808D+00 
-5 . 3267D+00 

8 . 0912D+01 
-4.13720+02 
-6 . 2475D+01 
-8.84200+00 

7.39920-01 
-2.4749D+01 
2.59050+01 
6 . 9321D-01 

7.2793D -02 
-1. 97820+00 
4.2967D+C0 
-3.1452D-01 

(b) 

Fixed throttle 






COLUMNS 
0. OOOOD+OO 
-4.5737D-01 
-2.7317D-02 
-1 . 3001D-02 

1 THRU 6 
0. OOOOD+OO 
*2 . 480 ID* 04 
-2.0654D+03 
-8.83480+02 

0. OOOOD+OO 
-7.28260+02 
-4.26320+00 
-1 . 3085D+01 

0. OOOOD+OO 
-1 . 21470+03 
8.01660+01 
1 . 1110D+01 

0. OOOOD+OO 
4 . 2798D+02 
4 . 9573D+02 
1.9376D+01 

0. OOOOD+OO 
-2.64490+02 
-2 . 0063D+01 
7 .46670+00 


COLUMNS 
0 . OOOOD+OO 
2 . 0012D+02 
2.5032D+02 
1. 2102D+01 

7 THRU 12 
0. OOOOD+OO 
-3 . 0735D+02 
7 , 0155D+00 
-1 . 0560D+01 

0. OOOOD+OO 
1 .54340+02 
-8.5866D+00 
-6 .15670-01 

0. OOOOD+OO 
-4.00260+03 
-3 . 3909D+02 
-1.4738D+02 

0. OOOOD+OO 
1.99500+02 
2 . 8377D+02 
1.40190+01 

0. OOOOD+OO 
3.6184D+01 
4.8903D+01 
2.1601D+00 


E-9 


Flight Condition: h*A0k ft, n*A 


(a) Throttle free 


COLUMN’S 
4.1124D-04 
■2.3227D-02 
4.3573D-03 
3 . 9995D-03 


1 THRU 6 
5 . 0036D+01 
-1.1466D+03 
“3 . 1217D+02 
-4 . 5419D+01 


1 . 5631D+00 
-1.68780+02 
1 . 0447D+00 
4 . 3753D+ 0 0 


-3 . 1841D+00 1.3668D+0Q 

-5 . 7583D+01 -8 . 95600+00 

1 .4759D+02 2.6414D+01 

7. 6267D+01 -1 . 0108D+01 


COLUMNS 7 THRU 12 


3.5106D-01 
-3. 3450D+01 
1 . 056 3D* 01 
- 3 . 8184D+00 


-4.7050D-02 
-1 .5111D+01 
6 . 9894D-01 
2.8959D-02 


2 . 00380-01 
-8 . 99720+00 
-4 . 0463D+00 
-6.4589D+00 


1 .5471D+01 
-3 .63140+02 
-1.0337D+02 
-1 . 7890D+01 


5.5993D-02 
1 .6913D+01 
2.85B4D+01 
1.54740*00 


(b) Fixed throttle 


COLUMNS 1 THRU 6 

0 . OOOOD+OO 0.00000*00 

-7.4833D-02 -3.1628D+03 

-3.1778D-03 -2.5123D+02 

2 .46170-03 -5.04060*01 

COLUMNS 7 THRU 12 

0 . 0O0CD*00 0 . O000D+00 

1.71930*01 -1.27550*02 

1.4775D*02 6 . 0470D*00 

-6 . 25480*00 -3.7364D-01 


0 . OOOOD* 00 
-3 . 0725D+02 
-1.03730*00 
-7. 36540-01 


0 . 0000D*00 
9.99490*00 
-6 . 7478D*00 
-8. *31680*00 


0 . OOOOD* 00 
-2 . 2405D*02 
5 .6016D+01 
5 . 1439D*01 


0 . 0000D+00 
-1.0546D*03 
-8.0806D*01 
-1.8416D*01 


0 . 00 OOD*00 
5. 2112D*01 
2 . 921 3D* 02 
-1 . 3608D+01 


0 . 0000D*00 
1.39090*01 
1.71520*02 
-2 . 9529D*00 


Flight Condition: h*50k ft, m-l.B, n*1 


(a) Throttle free 


COLUMNS 
-1.2428D-02 
-2.5982D-02 
-1 .5234D-05 
2 .66810-07 

COLUMNS 
3.1504D-01 
4. 4714D-01 
1.92140*01 
-1.8412D-01 


1 THRU 6 
1 ,4946D*03 
-1.22150*01 
4.55410-03 
1 . 5175D-05 

7 THRU 12 
-1 . 1959D*01 
-1 . 2251D*01 
-6.8022D-02 
-3.1292D-04 


-1.7242D*02 

-1.79270*02 

-1.6558D-01 

2.4938D-03 


-7.9929D-02 

8.1814D-02 

4,8313D*00 

-9.09480*00 


1.86570*04 
-1.7245D*Q1 
6 . 9925D-02 
-1. 15750-03 


-8 . 0056D*03 
-1 . 2989D*02 
-1.4003D-01 
2 . 1152D-03 


6.68810-01 
8.98050-01 
4.42740*01 
- 2 . 1426D+ 0 0 


3.8123D-01 
5. 18170-01 
2.41480*01 
-1.05170*00 


(b) Fixed throttle 


COLUMNS 1 THRU 6 


C.OGOOD*00 
-3. 9688D-01 
-8.5073D-04 
-6 .39870-06 

COLUMNS 
0 . OOOOD+OO 
1.4115D*02 
4.89440*01 
1.4580D*00 


0 . OOOOD+OO 
-3.8011D*04 
-8 . 1553D+01 
-6.1356D-01 

7 THRU 12 
0 . 0000D*00 
-1.65710*01 
-1.4756D-01 
-1.1025D-03 


0 . 00000*00 
-1.70610*02 
-3.4003D-01 
-2.5685D-03 


0 . OOOOD+OO 
-9.42180+00 
1.4088D+01 
-1.1752D+01 


0. OOOOD+OO 
-3.6535D+03 
-7.8179D+00 
-5.8789D-02 


0. OOOOD+OO 
0. OOOOD+OO 
0.00000+00 
0 . 00000+00 


0. OOOOD+OO 
2.5598D+02 
1 . 022 ID* 02 
1.4993D+00 


0. OOOOD+OO 
1.4298D+0 2 
5.6839D+01 
7.7799D-01 


5 . 9777D-01 
•1 . 7275D+01 
2.7939D+Q1 
2.56540+01 


1.5265D-02 
1 . 3191D+01 
5 . 1314D+00 
2.1119D-01 


0.00000+00 
-4 . 3588D+01 
-4.48070+01 
1.6665D+01 


0. OOOOD+OO 
2 . 5306D+CO 
2 . 9610D+C1 
-5 . 2277D-01 


2.9641D-01 
1 . 0008D-01 
-4.7749D+01 
1 . 29600+01 


6 .62080-02 
8.99770-02 
4 . 0977D+00 
-1.71690-01 


0. OOOOD+OO 
2.49610+01 
-5.77790+01 
1.0648D+C1 


0. OOOOD+OO 
2.4462D+C1 
9.7199D+00 
1. 3455D-01 


E-10 


Flight Condition: h=50k ft, m=1.8, n=2 


(a) Throttle free 


COLUMNS 1 THRU 6 



3 . 5959D- 03 
-2.8089D-02 
-2.7227D-03 
-6.3217D-04 

4 . 6794D+02 
-1 . 8184D+03 
-8 . 1445D+02 
-4 . 5145D+02 

3.5335D+00 
-1. 37080+02 
8 . 1206D+0 0 
7. 55180+00 

-1 . 1105D+01 
-1.4276D+02 
1.3217D+02 
1 . 0742D+02 

1.86490+00 
3 . 3696D+01 
2.20620+01 
-9.60430+00 

-8.0395D-01 

-3.50200+01 

-2.05620+01 

2.2742D+01 


COLUMNS ' 
4 . 6028D-01 
-1.6152D+01 
7.47510+00 
-4 , 0962D+00 

7 THRU 12 
-3.5162D-01 
-6.4189 D+00 
1 . 2381D+00 
6 . 7878D-01 

2.6473D-01 
-2.1512D+00 
-1 .74570+00 
-7.0306D+00 

7 . 7890D+01 
-3 . 1546D+02 
-1.4860D+02 
-8 . 6670D+01 

1.0763D-01 
3 . 0005D+01 
1.49240+01 
-1.4925D+00 

-1 . 1116D- 01 
7. 06400-00 
2.31070-00 
-7.20270-01 

(b) 

Fixed throttle 






COLUMNS 
0 . OOOOD+OO 
-6 . 0982D-02 
-3.6150D-03 
-1 .4051D-03 

1 THRU 6 
0 . OOOOD+O 0 
-4 . 1785D+03 
-5 . 4820D+02 
-2 .85590+02 

0. OOOOD+OO 
-2.5111D+02 
-9.6721D-01 
4.38100-01 

0. OOOOD+OO 
-2.82380+02 
3 . 8862D+ 01 
3 . 2729D+01 

0 .OOOOD+OO 
-6 . 9991D+ 01 
4 . 1818D+01 
-7.6396D+00 

0 . OOOOO+OC 
-7.16990+01 
-4 . 1017D+C1 
6. 26830+ OC 


COLUMNS 
0 .OOOOD+OO 
-7 . 29260+01 
1.8956D+01 
1.5 019D+00 

7 THRU 12 
0 . 0000D+ 00 
-3 . 0622D+01 
6 . 3265D-01 
-9 . 5266D-02 

0.00000+00 
1 . 9139D+01 
2.6529D+00 
-1 . 8004D+00 

0. OOOOD+OO 
-7 . 2199D+02 
-9.4700D+01 
-5 . 0892D+01 

0. OOOOD+OO 
-3.7145D+01 
2.4460D+01 
-1.3886D-01 

0. OOOOD+OO 
-5. 97160+ C' 
4. 16710+ CC- 

-1.68120-0: 

Flight 

Condition: h=50k 

ft, m= 1 . 8 




(a) 

Throttle 

free 






COLUMNS 
-8.7236D-06 
-1 .47060-02 
9 . 5740D- 03 
6.8283D-03 

1 THRU 6 
3.3564D+01 
-9.3343D+02 
-2 . 8197D+02 
-8.9545D+01 

1 . 3407D+00 
-1. 49340+02 
-1 . 1045D+00 
-2 . 3177D+00 

-2 . 8980D+00 
-3 . 1412D+01 
2 . 3251D+02 
1.4395D+02 

5.88410-01 
3.4538D+01 
1 . 3507D+01 
-1.5042D+01 

-2.1908D-01 
-2. 74940+ Cl 
-2 . 1212D+01 
3 . 3745D+01 


COLUMNS 
3 . 9985D-02 
-7.0423D+00 
5 . 1013D+00 
-4.87780+00 

7 THRU 12 
-6.0302D-02 
-1.6189D+00 
7 . 2874D-01 
2 .6307D-03 

5.0477D-02 

-5.6683D+00 

-3.3023D-01 

-7.9269D+00 

1 . 1111D+01 
-3.04290+02 
-1 . 0184D+02 
-3 . 6055D+01 

-1.18160-01 
5 . 0324D+01 
2.4240D+01 
8 . 3040D-01 

-8.0972D-03 
1 . 7827D+01 
4.4423D+0C 
4. 2546D-01 

(b) 

Fixed throttle 






COLUMNS 
0 . 00000+00 
-1 .77210-02 
1. 09310-03 
1.23980-03 

1 THRU 6 
0 . 0000D+00 
-9.5730D+02 
-8.2717D+01 
-4.0076D+01 

0. OOOOD+OO 
-1.7563D+02 
-2 . 1945D+00 
-3. 3685D+00 

0. OOOOD+OO 
-9 . 2676D+01 
4.2145D+01 
3.5458D+01 

0. OOOOD+OO 
5.93260+01 
4 .40640+01 
-7 . 1816D+00 

0. OOOOD+OO 
-2.86690+01 
-6 . 5929D+01 
7 . 5126D+00 


COLUMNS 

0.00000+00 

8.9326D+00 

1.99450+01 

-4.1092D-01 

7 THRU 12 
0 . OOOOD+OO 
-1.3331D+01 
4.7141D-01 
-4. 00340-01 

0. 0000D+00 
9 . 2995D+00 
6 . 5710D+ 00 
-1 . 8719D+00 

0.00000+00 
-3.4701D+02 
-2 . 9380D+01 
-1 . 5270D+01 

0. OOOOD+OO 
3 . 1186D+01 
2 . 7616D+01 
-2 . 1025D+00 

0. OOOOD+OO 
8 . 2989D+00 
4 . 7279D+00 
-1 .61540-01 


E— 1 1 



